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Abstract 


Microfluidic  geometries  allow  direct  observation  of  microscale  phe¬ 
nomena  while  conserving  liquid  volumes.  They  also  enable  modeling 
of  experimental  data  using  simplified  transport  equations  and  static 
force  balances.  This  is  possible  because  the  length  scales  of  these  ge¬ 
ometries  ensure  low  Re  conditions  approaching  the  Stokesian  limit, 
where  the  fluid  flow  profile  is  laminar,  viscous  forces  are  dominant 
and  inertial  forces  are  negligible.  This  work  presents  results  on  two 
transport  problems  in  microfluidic  geometries.  The  first  examines 
heterogeneous  binding  kinetics  in  a  microbead  array,  where  beads 
with  different  chemical  functionalities  are  sequentially  captured  in  a 
well  geometry  over  which  analyte  solution  is  flowed.  Finite  element 
simulations  identified  the  flow  rates  and  microbead  surface  receptor 
densities  at  which  the  binding  rate  approaches  the  kinetic  limit,  val¬ 
idating  the  results  for  the  prototype  NeutrAvidin-biotin  assay.  The 
second  part  of  this  work  discusses  the  dielectrophoretic  motion  of 
surfactant-stabilized  water  droplet  pairs  in  a  microchannel  as  they 
approach  and  coalesce  under  a  uniform  electric  held.  Experimental 
data  measuring  droplet-droplet  separation  distance  versus  time  were 
fitted  to  a  model  using  the  quasi-static  force  balance  between  the  at¬ 
tractive  electrostatic  force  and  the  resistive  hydrodynamic  force  with 


a  single  adjustable  parameter  representing  the  drag  force  coefficient 
between  each  droplet  and  the  adjacent  microchanncl  walls.  For  glass 
microchannels,  the  drag  force  coefficient  values  demonstrate  no-slip. 
However,  PDMS  microchannels  have  significantly  lower  coefficient  val¬ 
ues  corresponding  to  hydrodynamic  slip  lengths  of  1-2  /mi.  These 
large  slip  lengths  demonstrate  that  nanoporosity  plays  an  important 
role  in  the  hydrodynamics  of  PDMS  microchannels. 
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1 

Introduction 

1.1  Background 

One  of  the  main  themes  in  the  evolution  of  scientific  research  has  been  miniatur¬ 
ization.  From  early  advances  in  optics  that  allowed  direct  observation  of  micro¬ 
scopic  entities  such  as  cells  to  the  continued  progress  in  semiconductor  fabrication 
that  enables  increasingly  fast  (and  increasingly  portable)  electronics,  technolog¬ 
ical  innovation  is  often  synonymous  with  the  ideas  of  making  and  observing  on 
a  progressively  smaller  scale.  The  ramifications  of  this  miniaturization  are  man¬ 
ifold,  but  can  almost  always  find  application  at  much  larger  scales.  Investigation 
of  individual  cells  provides  insight  into  the  biology  of  living  organisms,  manu¬ 
facture  of  microscale  electronics  enables  communication  over  many  kilometers, 
and  the  synthesis  and  characterization  of  nanoparticles  improve  the  properties  of 
bulk  materials.  Fabrication  and  observation  of  microscopic  phenomena  are  in¬ 
valuable  tools,  but  in  many  instances  the  two  are  not  performed  simultaneously, 
precluding  the  opportunity  of  manipulating  the  microscopic  entities  in  real  time 
and  consequently  limiting  the  applicability  of  such  experiments  to  macroscale 
dynamic  processes.  To  make  a  microscale  observation  as  valid  as  possible  at  the 
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macroscale,  it  should  be  a  direct  miniaturization  of  the  macroscale  system,  in¬ 
cluding  all  of  its  components  and  capturing  all  of  its  dynamics. 

Such  a  miniaturized  system  can  be  called  a  lab-on-a-chip  (LoC),  implementing 
a  macroscale  process  in  a  microscopically  observable  environment  and  exploiting 
the  advantages  of  its  small  length  scale  to  reduce  the  time  and  increase  the  preci¬ 
sion  of  the  process.  The  LoC  concept  was  first  introduced  by  the  use  of  microelec¬ 
tromechanical  system  (MEMS)  components  in  electronics.  These  components, 
such  as  actuators,  gyroscopes  and  piezoelectric  elements,  are  fabricated  using 
the  same  photolithography  techniques  pioneered  by  the  semiconductor  industry 
to  manufacture  computer  processors  and  circuitboards.  Although  integration  of 
such  components  into  a  single  chip  has  become  ubiquitous  in  electronics,  the  LoC 
concept  has  been  popularized  by  the  use  of  photolithographically  patterned  ma¬ 
terials  to  act  as  templates  for  fluid  flow  geometries. 

Photolithography,  the  patterning  of  a  substrate  using  light,  has  evolved  sub¬ 
stantially  since  its  origins  several  decades  ago,  but  the  two  principal  components 
critical  to  its  success  as  a  pattern  transfer  technique  have  always  been  optics 
and  materials.  The  optical  path,  and  the  wavelength  of  the  light  used,  have  a 
dramatic  impact  on  the  maximum  achievable  fidelity  and  resolution.  To  achieve 
feature  sizes  on  the  order  of  nanometers,  currently  the  resolution  limit  for  semi¬ 
conductor  fabrication,  a  light  source  with  a  comparable  wavelength,  such  as  an 
electron  beam,  or  a  high  intensity  discharge  lamp  (KrF  or  ArF)  coupled  with 
an  optical  stepper  must  be  used.  The  tremendous  cost  of  this  equipment  limits 
the  availability  of  such  capabilities  to  a  small  number  of  institutions.  However, 
micrometer  resolutions  can  be  achieved  by  UV  exposure  sources  without  the  need 
for  an  optical  stepper  or  other  expensive  equipment.  This  resolution  limit  may 
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be  inadequate  for  applications  where  the  patterned  feature  density  demands  the 
highest  possible  resolution,  such  as  microprocessors  or  MEMS,  but  is  acceptable 
for  the  fabrication  of  fluidic  channels,  particularly  when  coupled  with  the  com¬ 
parable  resolution  of  bright  held  and  fluorescent  microscopy  methods  commonly 
used  to  observe  such  systems. 

The  materials  selected  for  the  fabrication  of  MEMS  and  semiconductor  com¬ 
ponents  are  dictated  by  precise  consideration  of  their  physical  properties.  In 
contrast,  the  materials  used  to  construct  fluidic  channels  are  in  most  applications 
not  critical  to  the  process,  provided  they  are  chemically  inert  to  the  huid(s)  with 
which  they  come  in  contact,  and  are  therefore  selected  primarily  by  considera¬ 
tions  of  cost  and  ease  of  fabrication.  To  fabricate  fluidic  geometries  from  glass  or 
crystalline  materials,  an  ablative  technique  such  as  etching  is  required.  A  number 
of  etching  methods  exist,  but  they  are  either  expensive  (reactive  ion  etching)  or 
time-consuming  and  dangerous  (HF)  to  the  point  that  they  are  of  limited  utility 
in  the  fabrication  of  devices  that  may  only  be  used  once.  For  this  reason,  pho¬ 
tolithographic  fabrication  of  fluid  channels  was  not  common  until  the  late  1990s, 
when  polymeric  materials  were  used  in  the  development  of  simple,  inexpensive 
fabrication  techniques  known  as  soft  lithography. 

Soft  lithography  was  pioneered  by  George  Whitesides  in  the  late  1990s  ((2)) 
and  it  dramatically  reduced  the  cost  and  difficulty  of  fabricating  fluidic  channels. 
Rather  than  use  the  patterned  photoresist  him  as  a  protective  layer  to  prevent 
regions  of  the  underlying  substrate  from  being  etched  or  coated,  soft  lithography 
uses  the  developed  him  as  a  mold  over  which  liquid  monomer  is  poured  and  poly¬ 
merized.  The  resulting  solid  is  then  simply  peeled  oh  the  mold,  taking  advantage 
of  the  polymer’s  large  elasticity  (and,  if  necessary,  chemical  release  agents  applied 
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to  the  photoresist  surface  prior  to  the  application  of  the  polymer).  The  photore¬ 
sist  him.  typically  consisting  of  a  crosslinked  epoxy,  is  extremely  durable  and  can 
be  repeatedly  reused.  Soft  lithography  therefore  enables  researchers  to  rapidly 
fabricate  large  numbers  of  devices  from  a  single  patterned  photoresist  mold  at 
ambient  conditions. 

The  most  commonly  used  polymer  in  soft  lithography  is  polydimethylsiloxane 
(PDMS).  The  chemistry  of  the  dimethylsiloxane  enables  it  to  be  oxidized  to  form 
silanol  groups,  which  can  react  with  other  silanol  groups  in  the  presence  of  trace 
amounts  of  water  to  form  siloxane  bonds.  This  silane  chemistry,  and  its  elasticity, 
enable  PDMS  to  be  covalently  bonded  to  itself,  as  well  as  to  glass  and  silicon, 
simply  by  bringing  the  surfaces  into  conformal  contact  subsequent  to  oxidation 
via  exposure  to  oxygen  plasma.  Cutting  access  ports  into  opposite  ends  of  the 
PDMS  layer  prior  to  bonding  allows  fluid  to  be  pumped  into  the  fluidic  chan¬ 
nel.  The  transparency  of  the  PDMS  to  optical  light  allows  direct  observation 
via  brightfield  microscopy,  and  the  low  fabrication  cost  allows  the  devices  to  be 
discarded  after  a  single  use. 

Fluid-based  LoCs  have  the  same  advantage  as  MEMS  in  terms  of  integrat¬ 
ing  multiple  components  (or  simply  multiple  channels).  Additionally,  the  small 
dimensions  of  the  device  require  only  a  microliter  scale  volume,  facilitating  the 
conservation  of  potentially  expensive  analytes,  ffowever,  the  single  most  impor¬ 
tant  advantage  of  these  fluid-based  LoCs  is  the  ability  to  exploit  the  inherent 
properties  of  transport  phenomena  on  small  length  scales  to  rapidly  reproduce 
experimental  conditions  that  are  easily  characterized  and  measured.  The  study 
of  these  phenomena  in  such  devices  is  known  as  microfluidics. 

The  most  fundamental  transport  process  occurring  in  the  microfluidic  channel 
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is  the  fluid  flow.  The  ability  to  precisely  control  this  flow  and  characterize  its 
hydrodynamics  is  central  to  accurate  microfluidic  experimentation.  Fortunately, 
this  is  easily  achieved  because  the  small  length  scale  ensures  microchannel  flows 
are  well  within  the  laminar  flow  regime.  Consequently,  fluid  streamlines  in  the 
channel  do  not  cross,  but  rather  different  fluid  layerss  flow  smoothly  past  one 
another  in  the  classic  conceptual  sense.  This  allows  the  generation  of  stable, 
smooth  interfaces  between  different  fluids  flowing  in  the  same  channel,  miscible 
or  immiscible.  It  also  has  important  consequences  with  respect  to  mass  trans¬ 
fer,  because  it  makes  mixing  in  microchannels  extremely  difficult.  Not  only  are 
microchannel  flows  laminar,  many  approach  the  Stokesian  limit,  in  which  inertia 
is  zero.  This  phenomenon  allows  the  researcher  to  actuate  the  fluid  flow  with  a 
precision  not  possible  in  macroscale  flows  like  a  faucet  or  a  river  by  simply  equi¬ 
librating  inlet  and  outlet  pressures  across  the  microchannel  to  instantaneously 
stop  the  flow  (in  practice  this  is  difficult  due  to  limits  in  instrumentation).  In 
theory,  inertia  is  always  present  and  the  fluid  must  decelerate  to  stop,  but  the 
time  scale  for  the  deceleration  is  vanishingly  small.  No  strict  definition  exists  for 
the  Reynolds  number  below  which  inertia  can  be  neglected,  but  this  assumption 
greatly  simplifies  the  hydrodynamics  of  particles  translating  (and  rotating)  in  the 
microchannel. 

Neglecting  inertia  simplifies  microchannel  hydrodynamics  to  the  Stokes  equa¬ 
tion,  but  analytical  solutions  can  still  be  difficult  to  obtain  due  to  the  microchan¬ 
nel  geometry.  This  is  partially  due  to  the  binary  nature  of  a  single  photoresist 
patterning  step,  resulting  in  a  geometry  that  requires  a  two-dimensional  solution. 
Such  a  solution  exists  for  the  simple  case  of  an  open  rectangular  duct  channel, 
which  simplifies  to  the  well-known  one  dimensional  solution  when  the  aspect  ra- 
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tio  of  the  duct  cross-section  is  sufficiently  large.  Many  microchannels,  however, 
contain  contractions,  expansions,  obstacles  and  other  geometries  that  complicate 
the  hydrodynamics.  These  cases  require  a  numerical  solution  using  finite  element 
or  another  computational  method.  Once  the  solution  for  the  fluid  velocity  profile 
has  been  obtained,  it  can  be  used  to  model  heat  and  mass  transport. 

Microfluidics  has  found  application  in  many  areas,  but  is  most  commonly  as¬ 
sociated  with  biological  research.  Numerous  reasons  for  this  exist.  A  low  cost, 
single  use,  biologically  inert  platform  that  is  readily  mounted  for  observation  us¬ 
ing  a  bright  field  or  confocal  microscope  and  has  theoretical  volume  requirement 
in  the  microliters,  microchannels  have  been  used  to  incubate  living  cells,  perform 
DNA  conjugation  and  amplification,  and  perform  bioassays.  An  entire  literature 
has  developed  in  the  course  of  less  than  two  decades  describing  different  ways 
to  exploit  the  novelty  of  microfluidic  transport  properties  to  enhance  biological 
research.  While  many  of  the  techniques  presented  in  a  biological  context  could 
in  theory  be  applied  to  other  areas,  material  compatibility  between  polymeric 
microchannels  such  as  PDMS  and  many  liquids  is  less  than  ideal.  Specifically, 
strong  solvents  such  as  toluene  and  chloroform  interact  strongly  with  the  PDMS 
material.  However,  as  discussed  earlier,  fabrication  of  microfluidic  devices  from 
glass  or  crystalline  materials  is  possible,  and  can  be  used  to  study  chemical  reac¬ 
tions  and  separations  in  microfluidic  geometries;  such  microreactors  are  the  focus 
of  many  studies. 

Microchannels  produced  by  soft  lithography  may  be  susceptible  to  solvents, 
but  they  can  still  be  used  with  some  non-polar  liquids,  specifically  higher  viscosity 
oils.  These  oils  are  inert  and  therefore  of  limited  interest  by  themselves.  However, 
their  use  as  a  carrier  fluid  for  a  dispersed  phase  is  now  a  common  area  of  research. 
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The  dispersed  phase  could  be  a  suspension  of  solid  particles,  and  numerous  stud¬ 
ies  of  the  hydrodynamics  of  such  systems  exist.  More  commonly,  these  oils  serve 
as  carrier  fluids  for  liquid  droplets  which  are  generated  in  the  microchannel  by 
intersecting  microchannel  flows  of  the  dispersed  phase  and  the  carrier  fluid  (con¬ 
tinuous  phase).  This  technique,  called  flow  focusing,  capitalizes  on  the  laminar 
flow  regime  in  the  microchannel  to  form  highly  monodisperse  droplets.  Typically, 
surfactant  is  required  at  the  liquid-liquid  (or  in  some  studies,  liquid-gas)  interface 
to  sufficiently  lower  the  surface  tension  to  allow  droplet  formation.  This  interface 
is  not  the  focus  of  study  in  many  applications,  and  the  droplets  themselves  are 
simply  carriers  for  an  analyte  (or  cell).  However,  the  importance  of  interfacial 
dynamics  in  scientific  fields  from  the  clinical  to  the  industrial  makes  research 
focusing  on  the  droplet  interface  extremely  interesting.  In  particular,  investiga¬ 
tions  of  droplet  coalescence  to  predict  emulsion  stability  and  droplet  wetting  on 
solid  surfaces  to  engineer  self-cleaning  surfaces  have  wide  application  and  have 
attracted  considerable  attention  in  microfluidics.  This  is  because  the  aforemen¬ 
tioned  advantages  of  microfluidics  naturally  lend  themselves  to  the  observation, 
manipulation  and  measurement  of  dynamic  processes  occurring  at  fluid-solid  and 
fluid-fluid  interfaces. 

1.2  Scope  of  Research 

This  study  will  demonstrate  the  ability  of  microfluidic  devices  fabricated  using 
simple  soft  lithography  techniques  to  act  as  platforms  for  the  quantitative  study 
of  interfacial  dynamics  in  two  applications.  The  first  is  a  heterogeneous  kinetics 
study  of  protein  screening  assays.  This  study  employed  a  microfluidic  geometry 
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to  fluidically  capture  and  retain  an  array  of  functionalized  glass  microbeads  to 
serve  as  substrates  for  a  prototype  heterogeneous  bioassay.  It  demonstrated  the 
ability  to  capture  and  retain  beads  at  addressable  locations  in  an  array  of  microw- 
clls  patterned  into  bottom  wall  of  a  PDMS  microchanncl.  This  retention  not  only 
allows  observation  of  the  microbead  surfaces  during  the  bioassay,  but  enables  the 
deposition  of  multiple  bead  sets  in  sequence.  By  recording  the  locations  of  beads 
from  each  set  as  they  are  successively  deposited  in  the  well  array,  an  address¬ 
able  registry  of  bead  surface  functionality  can  be  created.  This  eliminates  the 
requirement  to  identify  the  microbeads  using  optical  barcodes.  Once  completed, 
the  microbead  array  is  exposed  to  an  aqueous  solution  of  analyte  that  is  flowed 
through  the  channel.  The  analyte  binds  to  the  probe  molecule  displayed  on  the 
microbead  surfaces  according  to  classic  Langmuir  kinetics,  and  the  rate  of  binding 
is  determined  by  measuring  the  change  in  fluorescent  signal  on  the  bead  surfaces 
over  time  using  a  fluorescent  microscope.  The  use  of  large  n  arrays  of  different 
bead  functionalities  effectively  performs  precise  measurements  of  multiple  assays 
simultaneously  using  a  single  analyte  sample,  a  technique  known  as  multiplexing. 

To  interrogate  potential  ligand-receptor  pairings  by  observing  which  microbead 
surface  functionalities  the  analyte  molecule  binds  to,  or  to  test  a  sample  for  the 
presence  of  a  particular  analyte  using  beads  with  a  known  conjugate  receptor, 
precise  measurement  of  the  binding  rate  is  not  critical.  Rather,  a  limit  of  detec¬ 
tion  in  terms  of  the  fluorescence  intensity  on  the  bead  surface  must  be  exceeded  to 
identify  a  positive  analyte-surface  conjugation.  However,  the  ability  to  measure 
binding  rate  data  theoretically  enables  the  calculation  of  the  kinetic  constant  for 
the  reaction.  In  the  kinetic  limit,  the  analyte  concentration  is  constant  and  the 
surface  concentration  of  the  bound  analyte  as  a  function  of  time  can  be  easily 
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calculated.  However,  real  experimental  conditions  can  only  approach  this  limit, 
and  if  the  rate  of  reaction  at  the  surface  is  greater  than  the  rate  of  diffusion  to 
the  surface,  depletion  of  the  analyte  species  from  solution  will  occur.  The  result 
is  a  surface  reaction  in  which  the  analyte  concentration  in  the  fluid  layer  at  the 
microbead  surface  is  less  than  in  the  bulk  fluid,  and  fitting  such  data  to  calculate 
the  kinetic  constant  by  assuming  the  kinetic  limit  is  valid  would  produce  incorrect 
values.  Instead,  numerical  simulations  can  be  used  to  solve  the  fluid  velocity  field 
around  the  microbead  and  its  retaining  well.  This  time-independent  velocity  field 
can  then  be  coupled  to  the  unsteady  state  convective- diffusion  equation  to  solve 
the  mass  transfer  problem.  By  solving  the  concentration  profile  in  the  microflu¬ 
idic  geometry,  the  surface  concentration  of  the  bound  analyte  can  be  calculated 
as  a  function  of  time  and  the  result  compared  to  the  kinetic  limit.  By  varying 
the  flow  rate  of  the  incoming  analyte  solution,  as  well  as  the  reaction  rate  at 
the  microbead  surface,  experimental  conditions  at  which  the  binding  is  kineti- 
cally  limited  can  be  identified.  These  results  are  validated  for  the  well-known 
avi din-biotin  conjugation  using  glass  microbeads  conjugated  with  biotin  to  assay 
a  solution  containing  a  known  concentration  of  NeutrAvidin  (a  proprietary  form 
of  the  protein  avidin)  that  has  been  fluorescently  labeled  with  Texas  Red  dye. 
Performing  the  same  assay  at  different  flow  rates  allows  measurement  of  the  flu¬ 
orescent  signal  on  the  microbead  surfaces  in  each  experiment  to  be  compared  to 
the  simulation  results. 

The  second  part  of  this  work  is  a  hydrodynamics  study  of  droplet  electroco¬ 
alescence.  PDMS  microchannels  containing  a  flow  focusing  orifice  are  used  to 
generate  monodisperse,  surfactant-stabilized  water  droplets  in  mineral  oil.  The 
flow  rates  of  the  two  fluids  through  the  orifice  determine  the  droplet  spacing, 
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which  is  adjusted  so  that  the  droplets  are  propelled  through  the  channel  in  a 
single  hie  train  with  each  droplet  separated  from  adjacent  droplets  by  distances 
of  several  droplet  radii.  The  application  of  a  high-strength,  uniform  electric  held 
polarizes  the  conducting  water  droplets  in  the  insulating  mineral  oil,  and  when 
the  droplet  train  is  aligned  parallel  to  the  direction  of  the  applied  held,  adjacent 
droplets  experience  an  attractive  dielectrophoretic  force.  If  the  separation  dis¬ 
tances  between  all  droplets  in  the  train  are  equal,  then  the  attractive  forces  cancel 
and  the  droplets  experience  no  net  motion  relative  to  one  another  (neglecting  con¬ 
sideration  of  the  hrst  and  last  droplets  in  the  train).  However,  perturbations  of 
the  fluid  how  and  changes  in  orientation  of  the  droplet  train  relative  to  the  held 
as  the  droplets  translate  through  different  sections  of  the  microchannel  result  in 
droplets  pairing  oh  and  dielectrophoretically  moving  towards  one  another. 

The  relative  motion  of  the  two  droplets  as  they  approach  and  coalesce  can  be 
decoupled  from  the  translational  motion  of  the  droplets  through  the  microchannel 
due  to  the  pressure-driven  how  of  the  mineral  oil,  a  consequence  of  the  Stoke¬ 
sian  limit  near  which  the  hydrodynamics  occur.  The  trajectory  of  the  droplets 
is  therefore  determined  by  a  simple  force  balance  between  the  attractive  electro¬ 
static  force  and  the  resistive  hydrodynamic  drag  force.  Both  are  functions  of  the 
droplet-droplet  separation  distance,  which  can  be  precisely  measured  as  a  func¬ 
tion  of  time  using  a  high-speed  camera.  The  hydrodynamic  drag  force  is  also  a 
function  of  the  droplet  velocity,  which  is  the  time  derivative  of  the  separation  dis¬ 
tance.  Integration  of  this  form  of  the  force  balance  produces  a  model  expression 
for  the  droplet-droplet  separation  distance  as  a  function  of  time  that  can  be  fitted 
to  the  experimental  data.  If  the  hydrodynamic  drag  force  on  the  droplet  is  only  a 
function  of  the  droplet-droplet  separation  distance,  the  model  has  no  adjustable 
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parameters.  This  would  be  true  for  two  spheres  approaching  in  a  liquid  with 
no  additional  surfaces,  provided  the  boundary  condition  at  the  sphere  surface  is 
known.  The  surfactant  immobilization  of  the  oil-water  interface  produces  a  no¬ 
slip  boundary  condition,  so  the  droplets  can  be  treated  as  hard  spheres  and  the 
droplet-droplet  hydrodynamics  are  fully  defined.  However,  the  droplets  in  the 
microchannel  have  diameters  on  the  order  of  the  channel  height,  and  the  larger 
density  of  water  relative  to  the  mineral  oil  creates  a  negative  buoyancy  resulting 
in  droplet-wall  separation  distances  of  less  than  1  /mi.  This  introduces  an  addi¬ 
tional  component  to  the  hydrodynamic  drag  force  on  the  droplet,  one  that  is  a 
function  of  a  separation  distance  that  cannot  be  directly  measured  in  this  exper¬ 
iment.  This  force  is  therefore  introduced  as  a  fitting  parameter  in  the  integrated 
force  balance. 

In  addition  to  being  dependent  on  the  boundary  condition  at  the  droplet 
surface,  the  droplet-wall  hydrodynamic  drag  force  is  a  function  of  the  boundary 
condition  at  the  microchannel  surface.  The  standard  no-slip  assumption  again 
results  in  a  well-defined  problem,  with  the  fitted  value  of  the  force  then  used  to  di¬ 
rectly  calculate  the  droplet-wall  separation  distance.  An  independent  estimation 
of  the  droplet-wall  separation  distance  can  be  made  by  calculating  the  trajectory 
of  the  negatively  buoyant  droplets  based  on  their  known  translational  velocity 
and  assumed  initial  position  in  the  center  of  the  channel  upon  formation.  This 
distance  can  be  compared  to  the  data  fitted  value  to  validate  the  no-slip  boundary 
condition  at  the  microchannel  surface.  If  the  surface  is  glass,  the  no-slip  assump¬ 
tion  is  valid.  However,  PDMS  microchannels  are  shown  to  have  droplet-wall  drag 
force  values  that  are  significantly  lower  than  expected  for  a  no-slip  wall.  This  is 
because  PDMS  is  not  an  ideal  solid,  but  rather  a  nanoporous  matrix  due  to  its 
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air  solubility.  Nanopores  on  the  microchannel  surface  that  are  open  to  the  fluid 
flow  enable  finite  fluid  velocities  at  the  surface.  These  velocities  are  defined  as 
the  product  of  the  velocity  gradient  at  the  surface  and  a  slip  length  determined 
by  the  chemical  and  physical  properties  of  the  surface  in  the  classical  definition  of 
hydrodynamic  slip.  Slip  lengths  are  usually  neglected  for  fluid  flows  over  smooth 
solid  surfaces  because  their  magnitude  is  on  the  order  of  nanometers  and  such 
distances  are  negligible  relative  to  continuum  flow  dimensions.  However,  porous 
surfaces  can  have  much  larger  slip  lengths  dne  to  the  ability  of  small  pores  to  re¬ 
tain  air  and  not  be  wetted  by  the  mineral  oil  from  the  microchannel.  The  surface 
therefore  maintains  a  Cassie-Baxter  state  that  is  only  partially  wetted  and  the 
slip  length  is  a  function  of  the  ratio  of  the  viscosities  of  the  two  phases.  This  work 
reports  measurements  of  slip  lengths  of  both  hydrophobically  functionalized  glass 
and  native  PDMS  surfaces.  As  expected,  the  glass  surface  has  a  slip  length  of 
zero  to  the  precision  of  the  measurement,  while  the  PDMS  surfaces  demonstrate 
slip  lengths  of  1-2  /mi  consistent  with  a  Cassie-Baxter  state. 
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Mass  Transfer  Study  of  a 
Prototype  Bioassay  in  a 
Spatially-Indexed  Microbead 
Well  Array 

2.1  Background 

The  field  of  microfluidics  has  received  more  attention  from  research  with  biological 
applications  than  from  perhaps  any  other  area.  The  rapid,  low-cost  fabrication  of 
biologically  inert  polymeric  devices  has  enabled  researchers  to  fabricate  LoCs  that 
perform  reactions,  separations  and  detection  assays  while  conserving  valuable  an¬ 
alyte.  The  optical  transparency  of  these  devices  has  allowed  for  unprecedented 
observation  of  these  processes  in  situ,  in  many  cases  permitting  quantification  of 
individual  biological  entities  such  as  cells,  DNA  or  proteins. 

One  of  the  most  important  areas  of  research  involving  biological  species  is  the 
screening  of  the  binding  interactions  of  proteins.  Proteins  control  biological  ac¬ 
tivity  by  selectively  binding  to  target  species  including  other  proteins,  peptides, 
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nucleic  acids  and  small  molecules.  The  high  specificity  of  these  binding  inter¬ 
actions  results  in  thousands  of  different  proteins  performing  unique  functions  to 
form  a  complex  network  that  can  be  mapped  to  give  greater  insight  into  biology 
at  the  cellular,  tissue  and  organism  levels.  This  mapping  of  protein  functions, 
known  as  proteomics,  allows  researchers  to  build  a  library  that  can  be  used  to 
rapidly  detect  abnormal  or  diseased  states.  Protein  screening  is  also  used  in  the 
development  and  validation  of  pharmaceuticals  for  targeted  drug  delivery.  This 
field  is  responsible  for  testing  candidate  drug  molecules  against  potential  receptor 
proteins  to  analyze  and  quantify  the  specificity  of  the  candidate  drug  in  selec¬ 
tively  binding  to  particular  protein  receptors  to  achieve  a  desired  biological  effect 
and  reduce  or  eliminate  harmful  side  effects.  Protein  receptors  also  play  a  vital 
role  in  the  intake  of  toxins  and  disease  vectors  to  tissues  and  cells,  and  therefore 
screening  assays  are  crucial  to  clinical  diagnostic  identification  of  disease  markers 
as  well  as  pathogen  detection  in  applications  including  environmental  surveillance 
and  food  monitoring. 

Due  to  the  large  number  of  potential  receptors  in  any  biological  system,  effec¬ 
tive  protein  screening  protocols  require  thousands  of  ligand-receptor  pairings  to 
be  interrogated.  Current  high-throughput  procedures  for  performing  these  tests 
typically  employ  one  of  two  approaches:  microarrays  or  bead  arrays. 

The  most  traditional  methodology  for  bioassays  is  the  96-well  plate  enzyme- 
linked  immunosorbent  assay  (ELISA)  (3,  4,  5).  This  platform  (also  termed  mi¬ 
crotiter  plates)  was  developed  in  the  1970s  and  is  still  commonly  used.  In  this 
approach,  a  plate  (typically  polystyrene)  is  used  as  a  substrate  to  attach  cap¬ 
ture  antibodies  to  the  surfaces  of  the  wells  arranged  in  a  rectangular  matrix.  An 
analyte  sample  is  dispensed  into  each  of  the  wells,  and  the  capture  antibodies 
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then  selectively  bind  with  antigens  present  in  the  sample  to  tether  them  to  the 
substrate.  The  amount  of  antigen  present  is  quantified  by  a  second  detection 
antibody  coupled  to  some  type  of  detectable  entity,  such  as  a  fluorophore,  which 
can  be  quantified  by  optical  microscopy.  To  enhance  the  capability  of  this  plat¬ 
form,  researchers  have  miniaturized  the  wells,  resulting  in  384,  1536  and  even 
9600  well  plates  using  etched  glass  or  silicon  substrates  (6).  While  this  method 
can  be  used  for  the  detection  of  one  or  several  analytes  in  many  bioassays,  the 
limitation  of  a  single  antibody  per  well  makes  conventional  well  plates  poorly 
suited  for  multiplexed  assays.  Recent  research  has  focused  on  approaches  to  cir¬ 
cumvent  this  limitation  by  using  robotic  spotting  or  lithographic  patterning  to 
produce  microarrays  of  different  antibodies  on  a  single  chip. 

Microarrays  offer  significant  improvements  over  well  plate  assays  (6,  7,  8,  9, 
10,  11,  12,  13,  14,  15,  16,  17,  18,  19).  By  presenting  multiple  antibodies  to  a  single 
analyte  sample,  they  effectively  reduce  the  required  sample  volume  and  decrease 
throughput  time.  Early  work  demonstrated  this  concept  by  simply  arraying  mul¬ 
tiple  receptors  in  each  well  plate  to  create  a  high-throughput  ELISA  (20).  More 
recently,  high-density  planar  chips  have  been  created  using  modified  inkjet  printer 
heads  and  pins.  While  inkjet  heads  and  pins  can  locally  deposit  protein  solutions 
to  create  arrays  on  the  order  of  microns,  the  smallest  dispensable  liquid  droplets 
limit  printable  array  resolution.  Lithographic  techniques  have  been  employed 
to  increase  array  densities  even  further.  Martin  et.  al.  employed  a  microcon¬ 
tact  printing  approach,  using  an  elastomer  stamp  patterned  using  soft  lithgraphy 
to  introduce  antibody  proteins  onto  an  aminosilane-functionalized  glass  surface 
(21).  Other  researchers  have  used  dip-pen  lithography  to  both  pattern  surfaces 
for  subsequent  deposition  of  proteins  and  deposit  proteins  directly  on  surfaces 
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(22,  23). 

In  spite  of  the  obvious  advantages  of  microarrays,  they  suffer  from  the  same 
drawbacks  inherent  to  2-D  platforms  such  as  well  plates.  These  devices  often 
produce  low  signal-to-noise  ratios  due  to  non-specific  adsorption  of  both  capture 
antibody  and  antigen  on  the  surface,  although  passivation  of  the  background 
with  polyethylene  glycol  can  remedy  this  issue.  Microarrays  have  lower  analyte 
sample  volume  requirements  than  well  plates,  but  significant  volume  may  still 
be  required  to  immerse  the  substrate  in  the  analyte  solution.  More  importantly, 
the  incubation  times  required  for  assays  on  planar  substrates  are  long  due  to  the 
heterogeneous  reaction  kinetics.  Finally,  although  microarrays  theoretically  offer 
a  platform  for  rapid  screening  of  many  antibodies,  the  large  array  densities  make 
indexing  the  identities  of  hundreds  or  thousands  of  potentially  unique  antibodies 
extremely  difficult.  An  alternative  approach  in  which  multiple  antibodies  can  be 
efficiently  indexed  for  a  high-throughput,  multiplexed  assay  is  therefore  desirable. 

Bead-based  assays  have  several  advantages  over  2-D  planar  assays.  By  using 
the  surfaces  of  small  particles  in  the  micron  size  range  as  the  platform  for  teth¬ 
ering  of  the  target  antigens,  the  high  surface-to-volume  ratios  of  the  particles 
can  be  exploited,  allowing  for  nearly  solution  phase  kinetics  to  reduce  incubation 
times  and  reducing  the  required  sample  volume.  Sets  of  beads  can  be  spectrally 
encoded  to  identify  a  particular  receptor  on  the  surface  of  each  bead  set,  allowing 
for  multiple  bead  sets  to  be  incubated  together  with  an  analyte  solution  to  pro¬ 
duce  a  multiplexed  assay.  To  identify  the  receptor  on  the  surface  of  a  particular 
bead  and  determine  to  which  antigen  it  has  conjugated,  the  bead  suspension  is 
flowed  through  a  narrow  channel  wherein  the  spectral  labels  of  individual  beads 
are  sequentially  excited  by  a  series  of  lasers  and  the  emission  spectra  recorded 
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and  deconvoluted.  This  technique  is  known  as  flow  cytometry. 

Flow  cytometry  is  a  commercially  available  technology  that  has  been  suc¬ 
cessfully  employed  to  perform  a  number  of  bioassays.  In  many  instances,  it  has 
demonstrated  superior  accuracy  and  sensitivity  compared  to  conventional  plate- 
based  protocols.  Early  experiments  utilized  polystyrene  microbeads  labeled  with 
a  single  fluorophore  to  perform  assays.  Stewart  and  Steinkamp  demonstrated 
that  such  a  system  could  be  used  as  a  standard  to  count  cells  in  blood  samples 
(24).  Later  work  has  interrogated  the  interactions  between  antibodies  and  im- 
munoglobin,  hepatitis  C  and  phospholipids  (25,  26,  27). 

The  main  advantage  of  flow  cytometry  is  its  multiplexing  ability  (28).  While 
early  research  did  not  exploit  this  capability,  the  demand  to  simultaneously  in¬ 
terrogate  multiple  analytes  will  become  more  critical  as  genomics  and  proteomics 
research  identify  more  potential  targets  for  study.  Currently,  commercially  avail¬ 
able  platforms  such  as  Luminex  xMAP  can  provide  up  to  500  spectrally  encoded 
bead  sets  for  multiplexed  assays  (29)  .  These  systems  have  seen  widespread  use 
over  the  past  decade,  and  recent  research  has  demonstrated  the  potential  to  in¬ 
crease  the  number  of  unique  labels,  or  barcodes,  still  further  while  improving 
upon  the  method  of  labeling  (30). 

However,  the  labeling  of  bead  sets  in  flow  cytometry  experiments  can  be  ex¬ 
pensive  and  cumbersome.  Although  commercial  platforms  are  available,  they  typ¬ 
ically  consist  of  polystyrene  beads  labeled  with  fluorescent  molecules.  Polystyrene 
beads  may  be  poorly  suited  for  specific  applications,  and  fluorescent  molecules 
undergo  photobleaching  over  time.  For  all  flow  cytometry  experiments,  calibra¬ 
tion  of  each  individual  bead  set  is  required;  this  takes  time  and  can  introduce 
significant  error  into  the  experiment.  Eliminating  the  requirement  for  spectral 
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labeling  of  bead  sets  would  therefore  dramatically  improve  the  efficiency  of  bead- 
based  bioassays. 

As  an  alternative  to  spectral  labeling,  bead  sets  can  be  deposited  onto  a  sur¬ 
face  to  create  a  spatially-indexed  array.  This  approach  combines  the  miniatur¬ 
ization  of  microarrays  with  the  enhanced  efficiency  of  bead-based  assays,  format 
which  incorporates  advantages  of  both  the  flat  microarray  and  the  flow  cytomet¬ 
ric/microbead  array  for  parallel,  high-throughput  screening.  (31,  32,  33,  34).  To 
accomplish  this,  the  beads  must  be  captured  and  sequestered  on  a  surface.  In 
general,  several  methods  exist  to  pattern  immobilized  microbeads  on  a  flat  sub¬ 
strate  surface.  Microbeads  can  be  deposited  by  gravity  from  a  solution  placed 
above  the  surface,  and  then  affixed  to  the  surface  by  using  electrostatic  interac¬ 
tions  (35,  36,  37,  38,  39),  covalent  bonding  (40,  41),  protein-ligand  binding  (42), 
an  adhesive  layer  (43),  or  by  transferring  a  pre-formed  array  of  beads  onto  an  ad¬ 
herent  surface  (44,  45).  However,  capturing  beads  by  gravity-settling  in  an  array 
of  wells  inscribed  on  a  surface  (a  well-plate)  presents  a  simpler  solution  because 
it  does  not  rely  on  bead/surface  interactions,  and,  by  properly  sizing  the  wells 
to  be  only  slightly  larger  than  the  microbead  diameter,  single  microbeads  can 
captured  at  the  array  (well)  location,  which  simplifies  the  tracking  and  correla¬ 
tion  of  screening  events.  Walt  and  collaborators  (46)  first  pioneered  the  trapping 
of  beads  with  surface  probes  in  wells  for  screening  applications  by  etching  wells 
into  the  tips  of  individual  fibers  of  a  fiber  optic  bundle  to  form  a  well-plate,  an 
approach  which  also  allowed  for  individual  readouts  of  fluorescently  labeled  bind¬ 
ing  events.  Incorporating  a  well-plate  filled  with  beads  into  a  microfluidic  cell 
can  be  undertaken  in  either  of  two  ways.  As  studied  by  Bau  et  al  (47,  48,  49), 
beads  are  first  trapped  in  the  wells  of  a  well-plate,  and  the  plate  is  then  incorpo- 
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rated  as  the  bottom  of  a  microfluidic  flow  cell.  Bau  et  demonstrate  that  the  flow 
through  the  cell  does  not  lift  the  beads  out  of  the  wells  as  long  as  the  flow  rate 
is  below  a  critical  value.  Bau  et  al  also  demonstrated  that  this  ex-situ  method 
of  bead  assembly,  because  it  allows  unhindered  access  to  the  array  during  the 
insertion  of  the  beads  in  the  wells,  can  be  used  to  position  the  beads  in  the  wells 
by  micromanipulation,  so  that  an  array  can  be  assembled  with  beads  displaying 
different  probes  with  the  probe  identity  at  each  array  position  known.  Bau  et 
al  also  showed  that  the  wells  can  be  loaded  by  random  deposition  from  solution, 
and  in  this  case,  to  display  beads  with  different  probes  in  the  array,  they  encoded 
the  beads.  Instead  of  ex-situ  assembly,  microbeads  can  also  be  assembled  directly 


Figure  2.1:  Idealized  schematic  of  the  assembly  of  a  microbead  array  by  the 
gravitational  settling  of  microbeads  into  wells  incorporated  as  the  bottom  of  a 
broad  channel  of  rectangular  cross  section  in  a  microfluidic  cell. 

into  an  array  in  a  microfluidic  cell  in  one  step  by  using  an  unfilled  well-plate  as 
the  cell  bottom  and  streaming  a  suspension  of  beads  through  the  cell  at  a  suf¬ 
ficiently  low  velocity  to  allow  individual  beads  to  be  captured  in  the  wells  due 
to  gravity  or  the  application  of  an  external  held  (see  Fig.  2.1  and  refs. (50,  51) 
who  also  demonstrated  the  use  of  the  array  for  a  binding  assay).  To  maximize 
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both  the  speed  and  the  efficiency  of  the  microbead  capture  in  this  format,  electric 
and  magnetic  fields  have  been  applied  to  charged  or  paramagnetic  beads  (respec¬ 
tively)  to  direct  the  beads  into  the  wells  (52,  53,  54).  Fluid  suction  has  also 
been  used  to  assist  in  the  bead  capture;  holes  placed  at  the  bottom  of  the  wells 
provide  a  liquid  path  from  the  channel  above  the  wells  to  drains  (see  McDevitt 
et  al  (55,  56,  57,  58,  59,  60,  61,  62,  63,  64)  and  Ketterson(65),  and  applying  a 
pressure  drop  across  these  holes  propels  the  beads  into  the  wells.  This  capture 
approach  also  increases  fluid  flow  around  the  beads  and  therefore  improves  mass 
transfer  of  analyte  during  the  subsequent  bioassay. 

This  study  uses  a  microfluidic  geometry  in  which  the  microbeads  are  intro¬ 
duced  into  the  device  in  a  fluid  suspension  and  captured  in  recessed  well  array  due 
to  gravity  (Fig.  2.1).  Our  objective  is  to  study,  both  theoretically  using  numerical 
simulation  and  experimentally  in  a  microfluidic  flow  cell  with  a  prototype  assay, 
the  mass  transfer  in  the  binding  of  a  protein  from  solution  to  ligand  molecules 
displayed  on  the  bead  surface.  The  results  of  the  analysis  can  be  used  to  construct 
guidelines  for  incubation  times  or  injection  volumes  and  flow  rates  to  ensure  a 
particular  level  of  binding  for  detection(66,  67,  68,  69),  or  to  define  kinetically 
controlled  regimes  in  studies  of  the  intrinsic  binding  kinetics  of  a  receptor-ligand 
pairs. 

In  the  standard  biosensor  geometry,  a  surface  patch  of  capture  probes  (length 
t  and  width  ts )  is  localized  in  a  rectangular  channel  of  width  w  and  height  h  with 
h  w  and  ts  ~  w.  The  convective  flow  of  the  target  analyte,  entering  the  flow 
cell  with  concentration  cQ  is  driven  by  either  a  pressure  gradient  (which  we  will 
consider  here)  or  electrokinetically  by  a  electric  potential  gradient.  For  h  <C  w 
end  effects  are  neglected  and  the  flow  can  be  considered  unidirectional  (in  the  y 
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direction)  and  only  a  function  of  z,  with  average  velocity  U ;  for  pressure  driven 
flow  vy(z)  =  —  <  1  —  4|  —  |  >.  The  transport  of  the  target  molecule  in  solu¬ 
tion  to  the  channel  wall  consists  of  diffusion  across  the  (parallel)  convective  flow 
streamlines,  and  kinetic  binding  of  the  target  to  the  probe  once  the  target  has 
arrived  to  the  sublayer  of  solution  immediately  adjoining  the  surface(68).  With 
h  <C  ts,  this  mass  transfer  is  principally  two  dimensional.  The  time  scale  for  a 
target  molecule  to  be  convected  along  the  patch  is  tc  =  £/U,  and  the  time  scale 

for  the  target  to  diffuse  across  the  channel  is  —  h2/T>  where  D  is  the  target 

tlh 

diffusion  coefficient.  The  ratio  of  these  scales,  tc/tr>  =  -  defines  a  Peclet  num- 

ber  (Pe  =  Uh/D ).  Typically,  h  ~  l()2/jm,  w  ~  103/xm  and  U  ~  102  —  l()4/jm2/s 
corresponding  to  flow  rates  Q  ~  1CT1  —  102yLt^ / min.  Target  proteins  or  smaller 
biomolecular  ligands  have  molecular  weights  of  order  103  —  104  and  corresponding 
diffusion  coefficients  of  ~  lC^/imP/s  so  that  Pe  is  large,  of  order  10  —  104.  If,  in 
addition  to  Pe  >  1,  the  sensor  patch  t  is  short  enough  such  that  t/h  <  Pe ,  then 
the  time  for  diffusion  across  the  channel  is  smaller  than  the  time  required  for  the 
target  to  move  over  the  patch  (tc  <  t®),  and  target  can  only  reach  the  surface 
through  a  boundary  layer  with  a  thickness,  which  increases  with  distance  down 
the  channel  but  is  always  smaller  than  h  with  the  target  concentration  outside 
of  the  boundary  layer  approximately  equal  to  the  inlet  concentration  cQ.  For 
Pe  >  1  and  t/h  <C  Pe,  the  boundary  layer  becomes  asymptotically  small  in  Pe 
everywhere  along  the  patch,  and  the  flow  in  the  boundary  layer  is  approximately 
linear  in  the  direction  normal  to  the  surface  (v^  «  6 (U/h)(z  —  h/ 2)).  The  bound¬ 
ary  layer  thickness  at  the  downstream  end  of  the  patch,  6,  can  be  estimated  as 
the  thickness  for  which  the  time  for  diffusion  across  this  thickness  to  the  patch 
is  equal  to  the  time  for  a  target  molecule  riding  at  a  distance  5  from  the  wall 
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5  f  i  jh  1  1/3 

to  reach  the  end  of  the  patch,  i.e.  —  ~  or  —  ~  <  — —  >  .  When,  for 

D  { U/h}  o  h  [  Pe  J 

Pe  >  1,  the  patch  size  is  large  enough  such  that  £/h  >  Pe,  then  the  diffusion 
times  across  the  channel  is  of  order  or  shorter  than  the  average  convective  passage 
time  along  the  patch  at  the  far  downstream  end,  and  the  boundary  layer  grows 
and  reaches  through  the  channel  cross  section,  depleting  the  bulk  concentration. 

At  the  patch  surface  the  diffusive  flux  is  equal  to  the  kinetic  rate  (per  unit 
area)  with  which  target  binds  to  the  surface  probe.  Kinetic  binding  is  a  bi- 
molecular  process  of  which  the  most  elementary  is  the  Langmuir  kinetic  scheme, 

<9r 

—  =  kacs  {Too  —  r}  —  kdT  where  T  is  the  surface  concentration  of  bound  target 
at 

and  Too  is  the  maximum  number  of  targets  which  can  bind  (per  unit  area),  cs  is 
the  sublayer  concentration  at  the  surface  and  ka  and  kd  are  the  association  and 
disassociation  rate  constants,  respectively.  The  equilibrium  surface  density  (Te9) 

r  /c  k  c 

is  -At  =  -  where  k  =  -DP.  During  the  binding  process,  the  sblayer  concen- 

Too  1  +  k  kd 

tration  initially  decreases  due  to  kinetic  binding,  but  at  later  times  increases  as 
the  surface  begins  to  saturate,  causing  the  kinetic  flux  to  decrease  and  the  diffu¬ 
sive  flux  to  repopulate  the  sublayer. 

For  Pe  >  1  and  £/h  Pe,  the  diffusive  flux  to  the  surface  through  the  bound¬ 
ary  layer  scales  as  — h - — ,  where  cs  is  the  sublayer  concentration;  equating 

o 

this  flux  to  the  the  maximum  kinetic  flux  defines  a  scale  for  the  sublayer  con¬ 


centration,  — 
Co 


1  +  Da 


l/h  \ 
Pe  j 


,  where  the  Damkohlcr  number  Da,  is  defined 


{£) 


1/3 


>  1 


k  r  h 

as  Da  =  - ^  .  In  the  limit  Pe  >  1  and  £/h  Pe,  when  Da 
(fast  binding  kinetics  relative  to  diffusion),  the  sublayer  concentration  tends  to 
zero.  This  mass  transfer  controlled  regime  has  been  studied  extensively  as  the 
entrance  region  problem  (70,  71),  with  analytical  expressions  for  T  as  a  function 
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of  t  and  the  distance  along  the  sensor  surface.  As  the  target  flux  to  the  surface  is 


controlled  solely  by  the  diffusive  mass  transfer,  the  characteristic  time  for  the  tar- 

Dc0 

get  to  bind  to  an  equilibrium  surface  density  (teq^)  is  given  by  teq^—^~  ~  reg  or 


teq,D  k 

Teq,V  =  7  „  7V7  ~  £ 


h2/ T> 


1  +  k  1  h 


—  [>  {Pe}  1  ,/3  where  £  =  and  Teq ^  denotes  a  nondi- 


mensional  completion  time  scaled  by  the  diffusion  time  across  the  channel  (fn). 
The  parameter  £  is  the  ratio  of  the  channel  height  h  to  the  adsorption  depth 
Too /c0,  the  distance  above  the  surface  which  contains  (per  unit  area)  enough 
target  to  saturate  the  surface.  For  t/h  >  Pe,  analytical  solutions  can  also  be  ob¬ 
tained  for  Da  — >  oo  for  the  fully-developed  concentration  profile(71).  When,  for 

’  t/h  \  1/3 


Pe  >  1  and  t/h  <C  Pe,  Dal - >  <  1  the  binding  kinetics  are  slow  relative 

l  pe  J 

to  diffusion,  and  the  sublayer  concentration  remains  at  the  inlet  concentration. 
(This  is  also  true  for  arbitrary  Pe  and  £ / h  if  Da  — >  0.)  The  process  is  only 
controlled  by  the  binding  kinetics,  and  — —  =  1  —  e~eDa^l+^T  where  the  charac- 

1  P 


eq 


k 


and  Teqtk  is  a 


teristic  kinetic  time  for  equilibrium  binding  is  Teg  &  ~  -  . 

eDa  (  1  +  k 

nondimensional  time  scaled  by  the  diffusion  time  txi  (see  Goldstein  et  al  who  have 
extended  this  analytical  solution  for  small  Da( 72,  73)).  For  intermediate  values 
of  Da  and  Pe  >  1,  analytical  solutions  can  be  obtained  for  T /Teq  <C  1  (71,  74). 
When  the  surface  concentration  is  not  negligible,  analytical  solutions  cannot  be 
obtained  because  of  the  nonlinearity  of  the  kinetic  equation.  For  Pe  3>  1  and 
£/h  <C  Pe  and  Da  of  order  one,  boundary  layer  (two  compartment)  models  in 
which  the  Langmuir  kinetic  equation  and  a  relation  equating  the  boundary  layer 
flux  to  the  net  kinetic  adsorption  are  integrated  either  numerically  in  time  for 
the  average  surface  concentration  on  the  patch(75,  76,  77,  78,  79,  80,  81,  82,  83) 
or  analytically(84).  Over  the  past  several  years,  numerical  solutions  by  finite  ele- 
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ment  or  finite  difference  solution  of  the  convective  diffusion  equation  for  the  target 
coupled  to  the  kinetic  exchange  at  the  patch  boundary  have  been  obtained  for  ar¬ 
bitrary  values  of  Pe,  Da  and  l  jh  to  obtain  the  surface  concentration  of  target  as  a 
function  of  time  and  distance  along  the  patch(66,  67,  68,  69,  71,  74,  84,  85,  86,  87), 
and  these  have  been  compared  with  the  two  compartment  model  solution  and  the 
results  of  binding  experiments  (see  for  example  (88,  89). 

The  mass  transfer  of  target  to  receptors  on  the  surface  of  a  bead  situated  in 
a  well  at  the  flow  channel  bottom  presents  a  more  complex  mass  transfer  than 
target  transport  to  a  patch  of  receptors  on  the  channel  surface,  and  has  not  been 
studied  in  the  detail  of  the  sensor  patch.  In  this  case,  target  streams  over  the 
top  half  of  the  bead  surface;  at  large  Pe  a  boundary  layer  does  develop,  but  the 
flow  is  attenuated  by  the  well  walls  and  the  target  is  not  streamed  as  directly 
over  the  probe  (bead)  surface  as  in  the  case  of  the  patch.  Our  object  in  this 
study  is  to  compute  the  surface  concentration  of  the  target  on  the  bed  surface  for 
arbitrary  Da  and  (large)  Pe  by  numerical  simulation,  and  to  assess  the  effects  of 
the  attenuated  flow  and  compare  to  the  transport  of  target  to  a  surface  patch  of 
probes  on  the  microchannel  wall  under  identical  conditions  (same  values  of  Da 
and  Pe).  The  avidin-biotin  binding  experiments  using  the  microfluidic  flow  cell 
microbead  array  will  also  be  undertaken  to  validate  the  regimes  drawn  by  the 
numerical  simulations. 

2.2  Transport  Simulations 

We  consider  first  the  mass  transfer  of  target  to  probes  on  the  surface  of  a  mi¬ 
crobead  situated  in  an  isolated,  circular  well  located  at  the  bottom  of  a  mi- 
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crofluidic  flow  channel  of  rectangular  cross  section.  The  values  for  the  geometric 
parameters  are  set  to  be  equivalent  to  the  experimental  design,  the  channel  height 
h  =  80  /mi,  the  well  depth  d—  50/irn  and  diameter  2 r  =  70  /mi.  The  microbead 
with  radius,  a— 42  /mi  is  positioned  along  the  axis  of  the  well,  and  is  recessed, 
as  in  the  experiments,  located  equidistantly  (±  4  /mi)  from  the  top  and  bottom 
of  the  well.  The  well  is  positioned  centrally  with  respect  to  the  side  walls  of 
the  channel,  with  the  well  axis  a  distance  w  from  the  walls.  A  (nondimensional- 
ized)  cartesian  coordinate  system  (lengths  scaled  by  h)  is  located  with  an  origin 
above  the  well  axis  at  the  center  of  the  channel,  with  y  along  the  flow  direction, 
5  perpendicular  to  the  bottom  wall  and  x  perpendicular  to  the  side  walls.  The 
computational  domain  is  closed  by  entrance  (upstream)  and  exit  (downstream) 
cross  sections  of  the  channel  located  a  distance  L  from  the  well  center. 

The  flow  of  the  analyte  stream  provides  the  convective  flow  setting  for  the 
mass  transfer  of  the  target,  and  is  described  first.  The  analyte  is  modeled  as  an 
incompressible,  Newtonian  fluid  with  the  density  p  and  viscosity  /j  of  water  (p= 1 
kg  m3  and  p  =  1CT1  kg  m_1sec_1)  independent  of  the  analyte  concentration.  The 
flow  through  the  channel  is  driven  by  a  pressure  gradient,  and  is  implemented  by 
assigning  a  uniform  velocity  U  in  the  y  direction  across  the  inlet,  and  a  zero  pres¬ 
sure  (relative  to  the  inlet)  across  the  exit.  The  flow  is  governed  by  the  continuity 
(mass  conservation)  and  Navier-Stokes  equations(90), 

V  •  v  =  0  (2.1) 

=  — Vp  +  V2v  (2.2) 

where  the  nondimensional  variables  V  and  V2  are  the  gradient  and  Laplacian 
operators  (scaled  by  h),  v  is  the  velocity  vector  (scaled  by  U),  p  is  pressure 
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(nondimensionalized  by  pU2),  and  t'  is  time  (scaled  by  convective  time  h/U)  and 


93  = 


pUh 


is  the  Reynolds  number.  For  the  experimental  flow  conditions,  typical 


p 

for  microfluidic  screening,  U  ~  102  -104  pm  sec-1,  the  flow  Reynolds  number  is 
of  order  10-2-  1,  and  therefore  the  flow  is  primarily  dominated  by  viscous  forces 
retaining  limited  inertial  effects.  The  continuity  and  Navier-Stokes  equations  are 
solved  in  cartesian  coordinates  with  the  inlet  and  outlet  conditions,  and  bound¬ 
ary  conditions  of  no  slip  on  the  interior  walls  of  the  channel  and  well  and  the 
bead  surface.  The  solution  is  obtained  numerically  using  finite  elements,  and 
time  marching,  implemented  with  the  COMSOL  Multiphysics  simulation  pack¬ 
age  (4.2),  using  both  triangular  and  quadrilateral  meshes.  In  the  absence  of  the 
well,  and  when  L  is  sufficiently  large,  the  (steady)  flow  at  the  origin  is  a  unidirec¬ 
tional  Poiseuillc  flow  through  a  rectangular  cross  section,  independent  of  y  and 
given  by  vy(x,  z).  The  distance  L  =3  X  103/im  is  taken  to  be  large  enough  so  that 
this  Poiseuille  flow  is  obtained  when  the  time  step  is  small  enough  and  the  mesh 
density  is  fine  enough,  and  this  provides  a  first  validation  of  the  flow  simulations. 
w  is  then  taken  large  enough  (w— 3  X  103/jm)  so  that  the  Poiseuille  flow  becomes 
independent  of  x  (at  the  origin),  so  that  the  side  walls  do  not  influence  the  flow 
at  the  well.  These  flow  simulations  (and  the  resulting  mass  transfer  simulations) 
are  therefore  in  the  absence  of  hydrodynamic  effects  associated  with  the  channel 
inlet,  outlet  or  side  walls.  When  the  well  is  unoccupied,  the  hydrodynamics 
is  an  open  cavity  flow,  as  shown  for  93  =  1,  in  Fig.  2.2(a)  for  the  y  component  of 
the  velocity  profile  (normalized  by  the  average  velocity  U)  as  a  function  of  5  for 
5;=0  and  y  =  0  (the  well  centerline),  and  x  =  0  and  y  =  28/80  corresponding  to 
a  location  inside  the  well  and  between  the  bead  and  the  well  wall.  Fig.  2.2(b)  is 
a  plot  of  the  magnitude  of  the  velocity  field  in  the  plane  x  —  0.  In  this  plane,  a 
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separatrix  streamline  clips  into  the  well  a  distance  5  of  ~  .25  on  the  well  axis  and 
separates  recirculating  flow  in  the  cavity  from  the  primarily  unidirectional  flow 
in  the  y  direction  in  the  channel  (note  the  change  in  sign  of  the  y  component  of 
velocity).  The  recirculation  consists  of  one  large  eddy,  as  would  be  expected  since 
the  aspect  ratio  of  the  well  d/(2r)  =  5/7  is  less  than  one,  and  consecutive,  op¬ 
positely  rotating  eddies  at  the  center  develop  for  deep,  rather  than  shallow  wells. 
The  y  component  of  the  velocity  on  the  well  axis,  which  is,  for  z  =  —1/2  approxi¬ 
mately  one  half  of  the  average  velocity,  decreases  exponentially  with  distance  —5 
into  the  well.  The  flow  pattern  in  the  presence  of  the  microbead,  also  for  93=1,  is 
also  shown  in  Figs.  2.2(b).  The  separatrix  streamline  in  the  x  =0  plane  is  forced 
upwards  by  the  bead,  and  a  circulation  develops  between  the  microbead  and  the 
well  wall,  although,  as  evidenced  by  the  magnitude  of  the  y  component  of  velocity 
at  the  off  axis  position  (y  =  28/80),  is  very  small  below  the  separatrix.  These 
flow  patterns  make  apparent  that  when  the  well  is  occupied  by  a  bead,  the  direct 
streamline  flow  in  the  channel  only  contacts  directly  the  microbead  surface  at  the 
top  of  the  microbead  where  the  separatrix  streamline  dips  along  the  microbead 
surface,  and  the  remainder  of  the  microbead  surface  is  contacted  by  a  very  slow 
recirculating  flow  which  separates  from  the  mainstream. 

Simulations  of  the  rate  at  which  targets  bind  to  the  probes  on  the  surface  of 
the  microbeads  in  the  wells  from  the  analyte  solution  streaming  over  the  beads 
is  obtained  by  solving  the  convective-diffusion  equation  (eq.  2.3)  for  the  mass 
conservation  of  the  analyte  in  solution  (in  Cartesian  coordinates).  This  is  done 
using  a  finite  element  numerical  simulation  with  forward  marching  in  time  that 
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was  implemented  with  the  commerical  software  package  COMSOL. 

r)n 

—  +  Pe  v  ■  Vc  =  V2c  (2.3) 

or 


In  the  above,  c  is  the  concentration  of  target  (non-dimensionalized  by  the  inlet 
concentration  c„),  v  is  the  steady  velocity  obtained  above,  and,  as  in  the  Intro¬ 
duction,  Pe  =  Uh/T>  and  time  is  scaled  by  the  diffusion  time  t®  =  h2 /D.  We 
assume  at  the  inlet  cross  section  that  the  concentration  of  target  is  uniform  (cG), 
and  the  distribution  has  relaxed  completely  at  the  exit  so  that  the  derivative  with 
respect  to  the  flow  direction  y  is  equal  to  zero.  Eq.  2.3  is  solved  with  these  condi¬ 
tions,  and  assuming  zero  flux  of  solute  on  the  interior  channel  and  well  surfaces, 
and  equating  the  diffusive  flux  to  the  kinetic  adsorption  at  the  microbead  surface. 


(n  •  Mbead  =  Da 

dt 
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dr 


=  £ 


1  +  k 


1  +  k 

(n  '  Vc)bead 


1  +  k 


(2.4) 
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where  n  is  the  outward  normal  to  the  microbead  surface,  cs  is  the  nondimensional 

cQh 

sublayer  concentration  and,  as  before,  £  =  (Too  is  the  maximum  surface  con- 

^  OO 

k  r  h 

centration  of  target),  Da  =  and  r  is  the  surface  concentration  scaled  by 

r  /c  A/  c 

the  equilibrium  concentration,  re(J,  where  — —  =  1  |  -  and  k  =  (ka  and  kd 


l  +  k 


kd 


are  the  adsorption  and  desorption  rate  constants).  The  surface  concentration  is 
a  function  of  the  position  on  the  bead  surface,  and  we  denote  by  f  the  average 
value  on  the  bead  surface.  The  mesh  and  time  step  are  refined  until  T(r)  is  in¬ 
dependent  of  the  mesh  density  and  the  time  step. 

In  nondimensional  form,  the  target  binding  T(r)/roo  is  a  function  of  the 
Damkohlcr  and  Peclet  numbers,  k  and  e.  In  the  prototype  assay  experiments  to 
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be  described  later,  the  binding  equilibria  is  nearly  irreversible  (k  1),  a  common 
characteristic  of  receptor-ligand  binding  interactions.  Therefore,  the  simulations 
are  performed  using  the  approximation  that  k  is  infinite.  The  parameter  e  scales 
the  overall  time  for  equilibration.  In  most  screening  applications  the  concentra¬ 
tion  of  the  target  is  low  enough  or  the  binding  capacity  large  enough  so  that  the 
adsorption  depth,  T0 0/c0  -  the  distance  above  the  surface  containing  enough  ma¬ 
terial  to  saturate  the  surface  per  unit  area  -  is  large  relative  to  the  channel  height 
h  so  that  eCl.  This  is  also  true  in  the  experiments,  and  we  set  £  =  0.016  in  the 
simulations  which  is  the  experimental  value.  We  first  examine  the  case  of  Pe  = 
10,  a  value  at  the  low  end  of  the  range  of  values  of  the  Peclet  number  in  microflu- 
dic  screening.  In  Fig.  2.3(a)  the  surface  concentration  of  targets  as  a  function  of 
time  (T(r))  for  Da  =  1,  10  and  102  for  binding  to  the  surface  of  a  microbead  in 
a  well  is  shown.  This  binding  rate  on  the  microbead  surface  is  compared  to  the 
binding  of  target  from  a  Poiseuille  flow  onto  a  circular  patch  of  probes  situated 
centrally  at  the  bottom  of  the  microchannel  wall  (z  =  —1/2),  and  with  a  radius 
equal  to  the  well  radius  r  and  with  the  binding  capacity  Too  and  kinetic  rate  ka 
identical  to  that  on  the  micobead  surface.  In  the  nondimensional  form  presented 
in  Fig.  2.3(a)  with  r  nondimensionalized  by  the  diffusion  time  scale,  increasing 
Da  corresponds  to  a  binding  experiment  in  which  the  kinetic  binding  rate  ka  is 
increased,  with  the  average  velocity  U,  concentration  c0,  binding  capactity 
and  diffusion  coefficient  T>  held  fixed.  For  both  the  circular  patch  and  the  mi¬ 
crobead  surface,  the  concentration  of  bound  target  increases  monotonically  with 
r,  and  as  Da  increases,  the  binding  rate  is  observed  to  increase.  The  binding  of 
target  to  the  surface  probes  is  a  transport  process  of  bulk  diffusion  to  the  surface 
followed  by  the  kinetic  step  of  target-probe  conjugation.  The  process  begins  as 
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target  in  the  sublayer  of  analyte  immediately  adjacent  to  the  probe  surface  binds 
to  the  surface,  depleting  the  sublayer  concentration  cs.  Depletion  continues  until 
the  surface  kinetic  rate  becomes  reduced  by  the  partial  saturation  of  the  surface, 
in  which  case  bulk  diffusion  repopulates  the  sublayer  until  the  sublayer  returns 
to  cD  (nondimensionally  to  one).  For  the  smallest  values  of  Da,  kinetic  exchange 
is  much  slower  than  bulk  diffusion,  and  the  sublayer  concentration  remains  rel¬ 
atively  uniform  around  the  microbead  or  above  the  patch,  eliminating  diffusion 
barriers.  In  this  limit,  the  average  surface  concentration  is  given  by  the  expo¬ 
nential  expression  r(r)/roo  =  1  —  e~eDaT .  This  ideal  kinetic  limit  represents  the 
fastest  rate  at  which  target  can  bind  to  the  surface,  and  this  limiting  envelope  is 
shown  in  Fig.  2.3(a).  For  Pe  =  10,  this  kinetic  limit  is  only  coincident  with  the 
numerical  simulations  for  the  patch  and  the  microbeads  for  Da  <  .1  (data  not 
shown).  As  Da  increases  to  values  of  one  and  larger,  the  kinetic  rate  increases 
relative  to  diffusion  and  this  reduces  the  concentration  of  target  in  the  sublayer  of 
analyte  immediately  adjacent  to  the  surface,  cs ,  to  values  less  than  c0,  creating  a 
diffusive  barrier  to  binding.  Since  the  sublayer  concentration  is  no  longer  equal  to 
the  farheld  bulk  concentration,  but  is  smaller,  the  numerically  simulated  mixed 
diffusive-kinetic  binding  rate  falls  below  the  ideal  kinetic  limit,  as  is  evident  for 
Da  =1,  10  and  102  in  Fig.  2.3(b).  For  increasing  Da,  the  sublayer  concentration 
decreases  and  this  has  two  consequences:  First  the  numerically  simulated  mixed 
binding  rate  increases  as  the  diffusive  flux  to  the  surface  is  greater  the  lower  the 
sublayer  concentration.  Second,  relative  to  the  ideal  kinetic  limit,  the  mixed  sim¬ 
ulated  binding  becomes  increasingly  slower  (Fig.  2.3(a))  since  the  kinetic  limit 
assumes  the  sublayer  concentration  is  equal  to  the  farficld  concentration.  Fig. 
2.3(a)  also  makes  clear  that,  because  target  binds  more  quickly  to  the  patch  in- 
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terface  than  to  the  microbead  surface,  the  diffusive  barrier  which  develops  around 
the  microbread  is  much  larger  than  the  diffusive  barrier  which  develops  over  the 
patch. 

The  diffusive  flux  to  the  surface  of  the  microbead  is  smaller  than  to  the  sur¬ 
face  of  the  patch  because  the  diffusive  transport  in  the  case  of  a  patch  is  entirely 
through  a  convective  boundary  layer,  while  diffusion  to  the  the  microbead  surface 
is  through  a  convective  boundary  layer  over  the  top  of  the  bead  exposed  to  the 
flow,  but  through  a  trapped,  slowly  recirculating  flow  which  surrounds  the  bottom 
part  of  the  bead.  The  bulk  concentration  fields  and  the  surface  concentrations 
provide  more  detail  and  insight  into  this  difference  in  mass  transfer  between  the 
two  geometries.  Consider  first  binding  to  the  circular  patch  of  probes;  the  con¬ 
centration  above  the  patch  in  the  plane  x  =  0  (perpendicular  to  the  microchannel 
wall,  along  the  flow  and  at  the  center  of  the  patch),  and  the  surface  concentration 
in  the  plane  5  =  -1/2  (the  channel  wall)  for  Da  =  10  and  for  three  nondimen- 
sional  times,  is  shown  in  Figure  2.3(b).  For  this  relatively  small  value  of  Pe 
(10),  the  characteristic  patch  size  in  the  flow  direction  ( t  ~  2r),  relative  to  the 
channel  cross  section  h  (2 r/h  =  7/8),  is  still  less  than  Pe,  and  as  discussed  in 
the  Introduction  for  binding  to  a  patch,  for  l/h  <  Pe,  a  boundary  layer  forms 
over  the  patch  and  extends  into  the  streaming  flow  but  does  not  extend  to  the 
opposite  end  of  the  channel,  as  is  clear  in  Figure  2.3(b).  The  concentration  in  the 
boundary  layer  above  the  patch  shows  the  initial  depiction  in  the  concentration 
next  to  the  surface  due  to  the  large  value  of  Da ,  followed  by  an  increase  in  the 
sublayer  concentration  as  the  surface  begins  to  saturate.  As  the  boundary  layer 
is  thinner  at  the  upstream  part  of  the  patch,  the  diffusive  flux  is  greater  at  the 
front  end  of  the  patch,  and  the  surface  concentration  increases  and  saturates  from 
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the  upstream  to  the  downstream  end  of  the  patch.  Figure  2.3(b). 

For  the  binding  of  the  target  to  the  microbead  surface,  the  concentration  field 
in  the  plane  x  =  0  and  the  surface  concentration  along  the  hemisphere  x  >  0 
(projected  onto  a  circle)  for  Da  =  10  and  for  same  three  times  as  depicted  for 
the  patch,  are  shown  in  Fig.  2.3(c).  Again,  the  characteristic  length  of  the  probe 
area  in  the  streamwise  direction  £  ~  2 a  divided  by  h  (£/h  =  1/2)  is  smaller  than 
Pe,  and  a  boundary  layer  forms  above  the  separatrix.  In  the  region  in  which  the 
separatrix  is  directly  attached  to  the  microbead  surface,  target  diffuses  directly 
through  the  convective  boundary  layer  to  the  surface,  and  the  diffusive  flux  is  the 
largest  and  the  binding  rate  to  the  surface  the  greatest.  This  resembles  the  trans¬ 
port  to  the  surface  of  the  patch.  At  the  upstream  and  downstream  parts  of  the 
well  where  the  separatrix  dips  into  the  well,  target  diffuses  through  the  convective 
boundary  layer  and  then  through  the  slowly  recirculating  liquid  surrounding  the 
lower  half  of  the  microbead  in  the  well  to  reach  the  bead  surface.  The  diffusion 
through  the  essentially  stagnant  liquid  reduces  the  diffusive  flux,  and  the  liquid 
in  the  well  quickly  becomes  depicted  of  target  for  this  relatively  large  value  of 
Da.  As  a  result,  while  the  binding  rate  at  the  top  part  of  the  microbead  surface 
increases  rapidly,  the  surface  concentration  along  the  lower  part  increases  much 
more  slowly  (cf.  the  projection  of  the  surface  concentration),  providing  an  overall 
reduction  in  the  average  rate  of  binding  (T(r))  compared  to  the  rate  of  binding 
for  a  patch.  The  depletion  of  target  in  the  stagnant  liquid  in  the  well  is  slowly 
replenished,  and  the  equilibration  takes  a  much  longer  time  relative  to  either  the 
patch  equilibration  or  the  ideal  kinetic  limit.  As  with  the  patch  geometry,  the 
binding  to  the  surface  is  asymmetric  with  respect  to  the  flow  direction.  The  top 
part  of  the  microbead,  in  contact  with  the  thinner  part  of  the  boundary  layer  and 
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having  the  larger  diffusive  flux,  has  a  greater  rate  of  binding  relative  to  the  down¬ 
stream  part  of  the  microbead,  in  contact  with  the  thicker  part  of  the  boundary 
layer  and  a  reduced  diffusive  flux.  As  Da  decreases  and  the  ideal  kinetic  limit  is 
approached,  the  effect  of  the  stagnant  layer  around  the  microbead  in  the  well  in 
decreasing  the  diffusive  flux  is  reduced,  and  the  microbead  and  patch  geometries 
show  similar  binding  rates.  For  Da  <  .1  (data  not  shown),  the  concentration  of 
target  in  the  stagnant  layer  is  approximately  the  farfield  target  concentration  due 
to  the  large  kinetic  barrier  (relative  to  diffusion),  and  the  binding  rate  becomes 
identical  to  the  kinetic  limit. 

When  the  Pe  number  is  increased  to  a  value  of  104  (Fig.  2.4),  the  transport 
picture  changes  significantly  for  both  the  patch  and  microbead  geometries.  As 
discussed  in  the  Introduction,  when,  for  large  Pe,  the  characteristic  streamwise 
length  of  the  probe  area,  t,  divided  by  the  channel  height  h  is  much  smaller 
than  Pe,  convective  boundary  layers  of  target  over  the  probe  surface  develop  and 
become  very  thin.  The  corresponding  diffusive  flux  of  target  through  the  layer 
becomes  much  larger  relative  to  order  one  Pe,  and  this  increases  the  binding 
rate  of  the  target  to  the  probe  surface.  In  addition,  because  of  the  enhanced 
diffusion  rate,  for  any  value  of  Da  (and  particularly  large  values),  the  sublayer 
concentration  of  target  adjoining  the  probe  surface  is  not  depleted  by  kinetic  ad¬ 
sorption  to  the  extent  that  it  is  when  Pe  =  0(1),  and  the  mixed  diffusive-kinetic 
binding  becomes  closer  to  the  ideal  kinetic  limit.  These  results  are  evident  in 
Fig.  2.4(a)  for  T(r)  which  shows  clearly  that,  for  Pe  =  104,  the  binding  rates  for 
both  the  patch  and  the  microbead  geometries  (Ijh  Pe)  are  much  faster  than 
for  Pe  =  10  (compare  Figure  2.3(a))  at  the  same  values  of  Da,  and  are  closer 
to  the  ideal  kinetic  limit,  and  Figure  2.4(b)  where  the  boundary  layers  are  much 


33 


2.  MASS  TRANSFER  STUDY  OF  A  PROTOTYPE  BIOASSAY  IN 
A  SPATIALLY-INDEXED  MICROBEAD  WELL  ARRAY 


thinner  and  depletion  less  evident  when  compared  with  Pe  =  10  and  Da  =  10 
(compare  Figure  2.3(b)).  Consider  in  particular  first  the  patch  geometry.  For  Da 
=  1,  the  mixed  diffusive-kinetic  binding  rate  is  on  the  kinetic  envelope.  This  is  in 
agreement  with  the  criteria  established  in  the  Introduction,  for  which  kinetically 


limited  transport  is  valid  for  streamwise  patch  lengths  t  satisfying  l/h  <C  Pe 

£/h  \  1/3 


'  \  1/3 

when  Dai  — -  >  <  1  (for  Da  =  1,  Da 


Pe  j 


Pe  J 


0.1).  For  the  larger  values 


( e/h} 1/3 

of  Da  in  Fig.  2.4(a),  Dai  >  >1,  the  kinetically  limited  criteria  is  not  satis¬ 

fied,  and  as  is  evident  in  the  figure,  the  mixed  diffusive-kinetic  patch  simulations 
are  below  the  ideal  kinetic  limit.  In  the  case  of  binding  to  the  microbead  surface, 
some  depletion  of  target  still  occurs  in  the  stagnant  liquid  surrounding  the  lower 
part  of  the  microbead,  due  to  the  slower  diffusive  transport  in  this  liquid.  The 
liquid  inside  the  well  at  the  upstream  side  appears  to  be  more  depleted  of  target 
compared  to  the  liquid  at  the  opposite  side,  and  correspondingly  the  binding  rate 
is  slower  on  the  lower  part  of  the  upstream  end  of  the  microbead  surface  relative 
to  the  bottom  part  of  the  downstream  end.  This  contrasts  with  the  case  of  Pe 
=  10  (Figs.  2.3(c))  in  which  depletion  and  binding  were  more  symmetrical.  One 
reason  for  this  asymmetry  may  be  due  to  the  fact  that  the  boundary  layer  along 
the  separatrix  at  the  upstream  side  of  the  well  is  very  thin  due  to  the  high  Pe  and 
follows  the  contour  of  the  separatrix  streamline  which  dips  into  and  out  of  tthe 
well  just  upstream  of  the  microbead.  As  a  result,  a  strong  (lateral,  y  directed) 
diffusive  flux  is  directed  to  the  microbead  surface  along  the  ascending  part  of  the 
streamline  (i.e.  the  part  that  moves  out  of  the  well  and  just  next  to  the  target¬ 
binding  bead  surface)  reducing  the  z  directed  flux  through  the  separatrix  required 
to  bring  target  to  the  lower  part  of  the  microbead  on  the  upstream  side.  When 
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the  Peclet  number  is  small  and  equal  to  10,  the  boundary  layer  is  much  thicker 
and  extends  well  above  the  separatrix  at  the  upstream  side  (Figure  2.3(c))  and 
lateral  diffusion  is  not  as  significant  and  transport  is  principally  in  the  downward 
z  direction.  In  any  case,  the  reduced  diffusive  flux  through  the  stagnant  layer 
surrounding  the  microbead  in  the  well  accounts  for  a  large  barrier,  and  kinetically 
limited  transport  is  only  observed  for  Da  of  approximately  0.5. 

The  important  conclusion  which  can  be  drawn  from  these  simulations  is  that 
the  diffusive  barriers  to  mass  transfer  of  target  to  the  surface  of  a  microbead  in  a 
well  at  the  bottom  of  a  microfluidic  channel  due  to  the  stagnant  layer  of  analyte 
in  the  well  can  significantly  inhibit  the  binding  rate  relative  to  the  transport  to  a 
patch  on  the  bottom  surface  of  the  channel.  As  such,  only  very  slow  binding  rates 
(low  Da)  or  large  throughputs  (high  Pe )  can  ensure  for  the  microbead  geometry 
linetically  limited  binding  conditions.  All  the  simulations  presented  here  are  for 
a  single  well. 

2.3  Experimental  Setup 

2.3.1  Device  design 

The  microfluidic  geometry  consists  of  an  open  duct  channel  with  lateral  dimen¬ 
sions  of  15  x  5  mm  and  a  height  of  100  /mi.  The  bottom  wall  of  the  channel  is 
a  flat  surface  populated  by  a  uniform  array  of  circular  wells  70  /mi  in  diameter 
and  50  /mi  deep.  The  lateral  pitch  of  the  array  is  250  /mi  and  sequential  rows  are 
offset  by  125  /mi.  The  channel  is  connected  to  four  ports.  Entrance  and  exit  ports 
are  located  at  opposite  ends  of  the  channel  to  control  microbead  deposition  rates, 
remove  excess  microbeads  from  the  surface,  and  introduce  analyte  solution  using 
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syringe  pumps.  Ports  located  on  the  sides  of  the  channel  are  used  to  introduce 
the  different  sets  of  microbeads,  with  each  set  of  beads  connected  to  a  dedicated 
port  to  avoid  reintroduction  of  beads  after  indexing. 

2.3.2  Device  fabrication 

The  microfluidic  devices  used  in  the  experiments  were  constructed  from  two  lay¬ 
ers  of  polydimethylsiloxane  (PDMS)  via  soft  lithography.  SU8  2050  negative 
tone  photoresist  (Microchem)  was  spin-coated  (Laurell)  on  3”  silicon  wafers  at 
the  spin  speed  specified  by  the  manufacturer  to  produce  uniform  films  of  100 
H m  thickness.  After  a  subsequent  soft-bake  step  to  uniformly  evaporate  solvent 
from  the  him  using  hot  plates,  the  unpolymerized  photoresist  was  lithographi¬ 
cally  patterned  using  contact  exposure  of  transparency  masks  (Pageworks)  with 
a  collimated  xenon  mercury  light  source  (OAI)  passed  through  a  360  nm  long 
pass  filter  (Omega  Optical).  The  use  of  the  filter  eliminates  UV  light  below  350 
nm  that  would  otherwise  overexpose  the  photoresist  due  to  its  higher  absorbance 
at  lower  wavelengths.  It  also  facilitates  accurate  calculation  of  the  the  exposure 
time  required  to  polymerize  the  photoresist.  Negative  tone  photoresists  utilize 
UV  light  to  initiate  a  free  radical  polymerization  reaction  in  the  exposed  regions 
of  the  him.  After  exposure,  the  him  is  heated  to  accelerate  and  finish  the  reac¬ 
tion.  The  photoresist  is  then  developed  to  remove  the  unexposed  regions  of  the 
photoresist  him  by  immersing  the  wafer  in  solvent.  The  developed  photoresist 
is  rinsed  and  dried  before  being  returned  to  the  hot  plate  at  an  elevated  tem¬ 
perature  to  ensure  permanent  adhesion  between  the  crosslinked  SU8  epoxy  and 
the  silicon  substrate.  The  resulting  pattern  contains  the  negative  relief  of  the 
microhuidic  geometry  and  is  robust  enough  to  serve  as  a  reusable  mold  for  the 
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casting,  polymerization  and  removal  of  the  PDMS  layers. 

The  PDMS  layers  (Sylgard  184)  are  fabricated  from  a  10:1  by  mass  ratio  of 
dimethylsiloxane  base  and  curing  agent.  The  two  components  are  mixed  and  de¬ 
gassed  (Thinky),  then  poured  over  the  SU8  molds  and  further  degassed  using  a 
vacuum  pump.  The  molds  were  incubated  at  65  °C  for  two  hours  to  polymerize 
the  PDMS  before  the  patterned  layers  were  cut  and  peeled  from  the  molds  using 
a  scalpel.  To  enable  the  microfluidic  channel  to  be  connected  to  external  fluid 
flows,  access  ports  were  added  to  the  upper  layer  containing  the  channel  using  a 
1.5  mm  biopsy  punch  (Harris  Uni-Core).  The  channel  was  then  sealed  by  bonding 
the  two  PDMS  layers  together.  This  was  accomplished  by  exposing  the  surfaces 
to  be  bonded  to  an  oxygen  plasma  for  30  s  (Harrick)  and  then  bringing  them  into 
conformal  contact  to  form  the  microchannel.  The  oxygen  plasma  reacts  with  the 
PDMS  to  form  silanol  groups  on  the  surfaces,  which  react  with  silanol  groups 
on  the  opposite  surface  to  form  covalent  siloxane  bonds  that  permanently  seal 
the  microchannel.  To  facilitate  observation  using  a  microscope  stage,  each  device 
was  mounted  on  a  standard  glass  microscope  slide  using  an  additional  plasma 
exposure  step  to  again  produce  siloxane  bonds  between  the  (glass  and  PDMS) 
surfaces. 

2.3.3  Microbead  functionalization 

Glass  microbeads  with  a  mean  diameter  of  42.5  /jrn  (Duke  Scientific)  were  used  in 
this  experiment.  They  were  selected  for  ease  of  functionalization  and  sized  so  that 
only  a  single  bead  could  occupy  each  well.  Glass  beads  have  the  additional  ad¬ 
vantage  of  a  large  density  difference  with  water  to  facilitate  gravity-based  capture 
in  the  wells  of  the  array.  The  beads  were  first  cleaned  in  an  aqueous  solution  of  4 
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%  NH4OH  and  4  %  H202  (w/w)  heated  to  70  °C  for  30  minutes.  The  beads  were 
then  washed  twice  with  deionized  water,  twice  with  ethanol,  centrifuged,  washed 
twice  with  chloroform  and  centrifuged.  They  were  then  suspended  in  a  5  111M  so¬ 
lution  of  aminopropyltrimethoxysilane  (APS)  in  chloroform  for  one  hour  to  graft 
amine  groups  to  the  surface.  To  remove  unbound  APS,  the  beads  were  washed 
and  sonicated  with  chloroform.  They  were  then  centrifuged,  washed  twice  with 
ethanol,  twice  with  water  and  twice  with  dimethylformamide  (DMF).  Depending 
on  the  desired  functionalization,  the  beads  were  then  suspended  in  a  1  mg/mL 
solution  of  NHS-PEG,  NffS-PEG-Biotin,  or  NHS-Flnorescein  for  one  hour.  To 
ensure  no  photobleaching  of  the  fluorescein  occnred,  the  vial  was  covered  with 
aluminum  foil.  The  beads  were  then  washed  three  times  with  deionized  water 
and  stored  in  a  refrigerator  at  4  °C  until  used. 

2.3.4  Microbead  capture  and  Spatial  Indexing  Without 
Encoding 

To  introduce  the  functionalized  beads  into  the  microfluidic  channel,  the  device  is 
submerged  in  deionized  water  under  a  vacuum  to  remove  air  bubbles.  Polyethy¬ 
lene  tubing  is  inserted  into  the  access  ports,  and  the  tubing  connecting  to  the 
entrance  port  is  attached  to  a  plastic  syringe  controlled  by  a  syringe  pump  (Har¬ 
vard  Apparatus).  The  bead  sets  are  suspended  in  water  by  magnetic  stirring  and 
drawn  into  250  and  500  /iL  glass  syringes  (Hamilton).  The  glass  syringe  is  then 
connected  to  the  loading  port  on  the  side  of  the  microchannel,  elevated  above 
the  device  and  positioned  with  the  needle  facing  down  so  that  the  beads  fall  out 
of  the  needle,  into  the  tubing  and  down  towards  the  microchannel  under  gravity. 
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Fluid  flow  in  the  channel  is  initiated  at  10  /iL / min  so  that  the  beads  entering  the 
channel  are  immediately  propelled  along  the  surface  of  the  well  array  and  begin 
depositing  into  the  wells.  Beads  that  are  not  captured  by  the  wells  accumulate  in 
the  exit  port  and  are  not  carried  into  the  exit  tubing  because  of  the  moderate  flow 
rate.  They  are  returned  to  the  channel  by  reversing  the  flow  direction  using  the 
syringe  pump  so  that  the  beads  are  directed  towards  the  entrance  port.  The  flow 
direction  is  switched  repeatedly  until  the  desired  well  occupancy  is  achieved.  The 
flow  rate  is  then  increased  to  100  /iL/ min  (in  the  forward  direction),  resulting  in 
the  propulsion  of  the  beads  out  of  the  exit  port  and  into  the  exit  tubing.  After 
these  excess  beads  have  been  removed,  the  partially  occupied  well  array  is  imaged 
to  index  the  locations  of  the  beads  from  that  bead  set  using  a  lOx  inverted  micro¬ 
scope  objective  (Nikon)  connected  to  a  CCD  camera  (Scion)  controlled  by  IrnageJ 
software  (NIH).  This  setup  produces  images  with  a  field  of  view  encompassing 
20  wells  (5x4).  The  second  bead  set  is  then  deposited  into  the  microchannel 
and  the  loading  process  is  repeated  until  the  all  of  the  empty  wells  are  occupied. 
Residual  beads  are  washed  off,  and  the  well  array  is  again  imaged  to  verify  the 
locations  of  beads  from  each  set.  The  imaging  is  done  in  epifluorescent  mode 
using  a  100  watt  mercury  source  and  a  B-2A  filter  block  (Nikon)  to  selectively 
excite  and  detect  the  fluorescein-labeled  and  Texas  Red-conjugated  beads. 

The  ability  of  the  well  geometry  to  capture  and  retain  multiple  bead  sets  to 
create  a  spatially-indexed  array  of  bead  functionality  without  microbead  encod¬ 
ing  is  demonstrated  in  Figure  2.5.  The  first  bead  set  introduced  into  the  device 
is  functionalized  with  fluorescein  as  discussed  earlier.  Multiple  beads  are  cap¬ 
tured  by  the  well  array,  but  vacant  wells  remain  to  allow  subsequent  capture  of 
the  second  bead  set.  Llncaptured  fluorescein  beads  remaining  in  the  device  are 
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successfully  removed  without  displacing  the  captured  beads.  The  locations  of 
the  captured  beads  are  then  recorded  as  shown  in  Figures  2.5  for  two  different 
locations  within  the  well  array.  The  second  bead  set,  functionalized  with  PEG,  is 
then  introduced  into  the  array.  These  beads  are  captured  in  the  remaining  wells 
to  complete  the  index  of  bead  functionality  according  to  well  location. 

2.3.5  Prototype  assay 

The  prototype  assay  conjugating  NeutrAvidin  protein  labeled  with  Texas  Red  flu- 
orophore  to  biotin-functionalized  beads  was  performed  under  different  Pe  values 
to  corroborate  the  results  of  the  finite  element  simulations.  The  binding  curves 
shown  in  Figure  2.6  represents  the  mean  normalized  intensity  values  of  the  indi¬ 
vidual  beads  in  well  arrays  for  experiments  performed  at  Pe  =  5600  and  Pe  =  56. 
These  data  points  were  compared  to  binding  curves  obtained  by  COMSOL  for 
the  specified  Pe  value  and  different  values  of  Da  to  determine  the  value  of  Da  at 
which  both  experimentally  measured  binding  curves  match  their  predicted  curves. 
The  adsorption  rate  constant  ka  corresponding  to  the  matching  value  of  Da  is 
ka  =  7  x  104M-1s-1. 

2.4  Conclusions 

We  have  demonstrated  a  spatially-indexed  microbead  array  via  sequential  deposi¬ 
tion  of  microbeads  into  a  recessed  well  geometry.  Finite  element  simulations  used 
to  identify  the  flow  conditions  and  microbead  receptor  density  required  for  the 
observed  ligand-receptor  conjugation  to  approach  the  kinetically-limited  reaction 
rate  are  identified  and  validated  using  a  prototype  assay  of  NeutrAvidin  binding 
to  biotin  on  the  microbead  surfaces. 
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(a) 


(b) 


Figure  2.2:  (a)  Velocity  in  the  y  direction  (normalized  by  the  average  velocity 
U)  as  a  function  of  z,  in  the  plane  x  =  0,  at  the  center  of  the  well  ( y  =  0) 
and  at  an  upstream  location  inside  the  well  and  between  the  bead  and  the  well 
wall  (y=28/80)  in  the  presence  and  absence  of  a  microbead,  (b)  Magnitude  of 
the  velocity  in  the  plane  x  =  0  inside  the  well  in  the  absence  and  presence  of  a 
microbead.  All  simulations  are  for  lHc=l. 
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Varied  Da,  Pe  =  10,  e  =  0.016 


(a) 


(c) 


Figure  2.3:  Target  binding  to  probes  on  a  circular  patch  on  a  microchannel 
wall  and  on  the  surface  of  a  microbead  in  the  well  for  Pe  =  10:  (a)  The  average 
nondimensional  surface  concentration  on  a  surface  patch  and  the  surface  of  the 
microbead,  F,  as  a  function  of  t  for  Da  =10,  102  and  103.  (b)-(c)  Target  concen¬ 
tration  boundary  layers  around,  and  the  spatial  distribution  along  either  a  surface 
patch  (b),  or  the  microbead  in  the  well  (c)  for  r  =  1,  15  and  30  and  Da  =  10.  For 
the  microbead,  the  concentration  boundary  layer  is  in  the  plane  x=0,  and  the  sur¬ 
face  concentation  is  the  projection  of  the  concentration  on  the  hemisphere  x  >  0. 
e  =  .016  and  k  — >  oo.  Simulations  at  additional  Pe  and  Da  values  are  presented 
in  Appendix  A. 
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Varied  Da,  Pe  -  10000,  e  =  0.016 


(a) 


(b) 
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Figure  2.4:  Target  binding  to  probes  on  a  circular  patch  on  a  microchannel  wall 
and  on  the  surface  of  a  microbead  in  the  well  for  Pe  =  104:  (a)  The  average 
nondimensional  surface  concentration  on  a  surface  patch  and  the  surface  of  the 
microbead,  T,  as  a  function  of  r  for  Da  =10,  102  and  103.  (b)  Target  concentration 
boundary  layers  around,  and  the  spatial  distribution  along,  either  a  surface  patch 
(b)  or  the  microbead  in  the  well  (c)  for  r  =  1,  5  and  10  and  Da  =  10.  For  the 
microbead,  the  concentration  boundary  layer  is  in  the  plane  x=0,  and  the  surface 
concentation  is  the  projection  of  the  concentration  on  the  hemisphere  x  >  0.  e  = 
.016  and  k  — y  oo. 
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Figure  2.5:  Bright  field  and  fluorescence  images  of  sequential  bead  array  in  two 
fields  of  view. 
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Figure  2.6:  Normalized  binding  curves  for  prototype  NeutrAvidin-biotin  assay 
compared  to  finite  element  simulation  results  at  equal  Pe.  =  4.2  x  10“9M . 
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3 

Hydrodynamic  Slip  at  an 
Oil/Polydimethylsiloxane  Surface 
From  Measurements  of  the 
Dielectrophoretic  Motion  of 
Water  Droplets 

3.1  Background 

Recently,  new  attention  (91,  92)  has  been  paid  to  the  possibility  of  hydrody¬ 
namic  slip  at  an  interface  between  a  stationary  solid  surface  and  a  “simple” 
(non-polymeric)  liquid  moving  over  the  surface.  The  consideration  of  boundary 
slip  began  with  the  continuum  level  formulation  of  the  Navier  slip  condition  for 
a  Newtonian  fluid,  which  equated  the  fluid  velocity  tangent  to  the  surface,  vs,  to 
the  boundary  tangential  stress,  rs,  by  the  slip  coefficient  A,  i.e.  vs  =  jTs  where  /j 
is  the  fluid  viscosity  and  A  has  units  of  length.  Since  this  formulation,  experimen¬ 
tal  studies  have  made  clear  that  the  “no-slip”  condition  of  A  =  0  is  sufficient  to 
accurately  model  most  macroscopic  flows  with  length  scales  in  the  range  of  mil- 
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limeters  to  meters,  with  the  exception  of  contact  line  motion  (93)  and  the  flow  of 
polymeric  (non-Newtonian)  fluids  (94).  However,  recent  molecular  dynamics  sim¬ 
ulations  on  atomically  smooth  surfaces  have  demonstrated  slip  on  the  molecular 
scale,  and  calculated  A  as  a  function  of  the  strength  of  the  liquid-solid  interaction 
(95,  96,  97).  For  strong  liquid-solid  interactions  which  characterize  complete  or 
strong  wetting  of  the  liquid  on  the  solid  surface,  slip  lengths  are  of  the  order  of 
only  a  few  molecular  diameters  (0(1  nm)),  while  relatively  weaker  interactions  of 
partially  wetting  fluids  have  slip  lengths  extending  tens  of  diameters  (0(10  nm)). 
Current  tools  for  measuring  A  include  particle  image  velocimetry  (PIV),  image 
velocimetry  enhanced  with  evanescent  near-field  illumination  at  the  surface,  and 
atomic  force  microscope  (AFM)  and  surface  force  apparatus  (SF)  measurements 
(98,  99,  100,  101,  102,  103,  104,  105,  106,  107,  108,  109).  Several  experiments  are 
consistent  with  the  MD  calculations.  For  smooth  surfaces  of  a  wetting  liquid,  e.g. 
water  flowing  along  a  hydrophilic  surface,  either  zero  slip  or  coefficients  less  than 
a  few  nanometers  are  recorded,  while  water  over  a  partially  wetting  hydropho¬ 
bic  surface  (e.g.  self-assembled  octadecyl  silane  (OTS)  monolayers)  obtain  slip 
lengths  of  tens  of  nanometers,  e.g.  (108,  110,  111,  112).  For  nonpolar  liquids 
wetting  smooth  hydrophobic  surfaces  (the  weakest  liquid/surface  interactions), 
slip  lengths  are  of  the  order  of  a  few  tens  of  nanometers  (113,  114,  115). 

“Giant”  slip  (A  ~  0(1  //m)  or  larger)  is  the  subject  of  great  interest  for  its 
application  to  reducing  surface  friction  in  micro  and  nanofluidic  channel  flows. 
Large  slip  can  be  achieved  when  low  friction  air  layers  are  situated  between  the  liq¬ 
uid  and  the  surface  a  circumstance  which  arises  when  a  population  of  nanoscopic 
gaseous  domains  adhere  to  a  surface,  or  air  becomes  trapped  in  a  micro  or  nano- 
textured  surface  that  is  not  wet  by  the  liquid  (super hydrophobicity)  (116).  AFM 
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Figure  3.1:  Measurement  of  microchannel  slip  at  an  oil/PDMS  surface  by  observ¬ 
ing  the  dielectrophoretic  merging  of  water  droplets  in  oil  moving  in  close  proximity 
to  the  PDMS  channel  wall. 


studies  have  provided  direct  evidence  of  nanoscopic  gaseous  domains,  primarily 
at  the  interface  between  water  and  a  hydrophobic  surface  (117).  AFM  measure¬ 
ments  have  also  verified  the  reduction  in  surface  friction  (118),  and  nanoscopic 
gas  domains  have  been  suggested  (119)  as  one  reason  why  some  measurements 
of  slip  at  the  water/hydrophobic  solid  surface  (120,  121)  obtain  one  micron  or 
larger  slip  lengths.  Hydrophobic,  textured  surfaces  filled  with  air  provide  a  more 
reproducible  method  for  generating  large  slip,  (122,  123,  124,  125),  and  theoret¬ 
ical  MD  and  continuum  studies  of  model  textures  (see  for  example  (126,  127)) 
demonstrate  that  large  slip  requires  the  solid  fraction  of  the  surface  to  be  a  few 
percent. 

Little  attention  has  been  paid  to  obtaining  giant  slip  when  nonpolar  liquids 
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slip  over  a  surface  (“oil  slippery  surfaces”),  although  significant  interest  is  devel¬ 
oping  dne  to  the  emergence  of  dropwise  microfluidic  platforms  which  are  based 
on  water  droplets  moving  in  a  continuous  oil  stream  (128).  Generating  gaseous 
domains  at  the  surface  is  the  key  to  large  slip,  and  one  method  for  generating 
a  significant  coverage  of  nanobubbles  on  surfaces  submerged  in  oil  is  to  nucleate 
them  spontaneously  on  oil  contact.  Spontaneous  nonequilibrium  formation  of 
nanobubbles  at  a  surface  occurs,  for  example,  when  an  air  saturated  liquid  (e.g. 
ethanol)  is  displaced  by  a  second,  miscible,  liquid  (e.g.  water),  also  saturated  with 
air,  but  with  a  lower  solubility  (117)  and  nanobubbles  nucleate  to  accommodate 
the  reduced  solubility  of  the  displacing  phase.  In  this  chapter  we  demonstrate 
that  a  nanoporous,  hydrophobic  polymer,  polydimethylsiloxane  (PDMS),  which 
has  a  significant  permeability  to  air  because  of  the  hydrophobicity  of  the  polymer 
and  a  relatively  large  free  volume  (129),  permits  oil  to  slip  over  its  surface  with  an 
order  one  micron  slip  length.  When  oil  contacts  the  surface  of  an  air-equilibrated 
PDMS  substrate,  air  is  released  as  nanobubbles  to  the  oil/PDMS  surface.  The 
driving  force  for  this  release  derives  from  the  fact  that  the  PDMS  is  also  perme¬ 
able  to  oil  because  of  its  hydrophobic  nanoporosity,  and  diffusion  of  oil  into  the 
nanopores  on  the  contact  of  oil  to  the  surface  displaces  the  air  to  the  interface. 

We  demonstrate  in  particular  oil  slip  at  the  inside  surface  of  a  microchannel 
formed  in  a  PDMS  monolith  fabricated  by  soft  lithography.  We  use  a  viscous 
mineral  oil  {non  ~  100 nH2o),  a  mix  of  high  and  low  molecular  weight  olefins, 
which  does  not  macroscopically  swell  the  PDMS  (130)  and  distort  the  surface 
from  the  atomically  smooth  topology  normally  evident  in  AFM  measurements 
but  can  still  displace  air  in  PDMS  by  the  solubilization  of  small  linear  alkanes 
in  the  oil.  The  microfluidic  arrangement  is  shown  in  Fig.  3.1.  (Additional  de- 
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tails  are  given  in  B.)  Using  optical  microscopy,  we  measure,  in  a  PDMS  channel  of 
rectangular  cross  section  with  height  h  (100  /an)  and  wide  width  w  (300  /an),  the 
edge-to-edge  separation  distance  s(t )  of  pairs  of  nearly  occluding  water  droplets 
(“A”  and  “B”)  which  are  entrained  in  a  mineral  oil  stream  and  are  driven  to¬ 
gether  by  a  dielectrophoretic  (DEP)  force  of  attraction,  Fdep •  This  force  is  due 
to  an  electric  field  £  applied  parallel  to  the  bottom  wall  of  the  channel  and  along 
the  flow  direction.  The  merging  droplet  pair  is  part  of  a  single  file  droplet  train, 
formed  upstream  by  flow  focusing  (131)  of  oil  and  water  streams  (flow  rates  equal 
to  0.4-1  /tty min  and  0.04-0.1  /if'/ ruin,  respectively),  through  an  orifice  50  fin l  in 
width.  The  focusing  forms  droplets  with  radius  a  equal  to  approximately  40  pm, 
which  are  separated  by  a  few  radii,  and  flow  at  the  average  stream  velocity  V  of 
approximately  250  -  500  /t m/sec.  The  electric  field  is  applied  to  switchback  flow 
lanes  at  the  downstream  end  of  the  chip  by  a  voltage  V  set  across  parallel  copper 
strip  electrodes  inserted  through  the  PDMS  (a  dielectric)  to  insure  a  uniform 
field  across  the  flow  channels.  Relative  to  the  aqueous  phase  which  is  de-ionized 
water,  the  oil  is  nonconducting  and  the  field  polarizes  each  of  the  droplets  of 
the  pair  into  dipoles.  In  the  flow  lanes  parallel  to  £,  the  polarized  drops  are 
aligned  with  £  and  attract  each  other,  creating  the  dielectrophoretic  force.  Prior 
to  the  application  of  the  field,  the  train  flow  is  observed,  and  when  a  pair  are 
observed  to  pair-off  to  a  relatively  close  separation  (less  than  one  radius)  due  to 
flow  disturbances,  the  held  is  applied  to  merge  the  pair,  and  a  video  recording 
is  made  with  a  high  speed  camera.  At  the  time  of  application  of  the  held,  the 
droplets,  heavier  than  the  oil  ( Ph2o  =  103  kg/m3  and  pou  =  8.75  xlO2  kg/m3) 
have  settled  to  a  separation  distance  d  from  the  bottom  wall  of  the  microchannel 
which  is  determined  by  their  settling  velocity  and  the  transit  time,  r,  from  their 
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formation  at  the  orifice  until  application  of  the  field  and  is  of  the  order  of  a  few 

hundred  nanometers.  At  these  distances  from  the  wall,  the  approach  velocity 
ds 

Va  —  Vb  =  —  due  to  the  DEP  force  is  affected  (as  we  show  by  numerical  so¬ 
lution  of  the  hydrodynamic  equations)  by  the  drag,  against  the  bottom  wall,  of 
the  intervening  oil  between  the  droplet  and  the  wall,  and  the  slip  on  this  wall. 
From  comparison  of  s(t )  with  numerical  solutions,  A  is  obtained  with  a  precision 
which  can  distinguish  a  micron  size  slip  length.  We  also  measure  the  slip  when 
the  bottom  surface  is  a  glass  slide,  functionalized  with  octadecyltrichlorosilane 
(OTS),  which  we  do  not  expect  to  nucleate  nanobubbles  at  the  surface  because 
of  low  air  permeability,  and  we  find  zero  slip  to  the  precision  of  the  measurement. 

One  illustrative  data  set  is  given  in  Fig.  3.2  (The  video  is  in  the  supplemental 
material).  An  edge  detection  routine  is  used  to  determine  the  droplet  perimeters 
from  which  the  size  of  the  merging  droplets  a  and  the  pair  separation  distance 
s(t)  are  computed.  The  separation  distance  is  shown  in  Fig.  3.2,  and  the  relative 
approach  velocity  is  approximately  40  /jrn/sec.  To  obtain  different  data  sets  of 
s(t)  corresponding  to  different  droplet  radii  or  droplet-wall  separations  d  (transit 
times  r),  the  oil  and  water  flow  rates  at  the  flow  focusing  orifice  are  changed,  or 
the  merging  at  different  switchback  lanes  is  observed. 

A  nonionic  surfactant  (Span  80,  sorbitan  monooleate)  is  dissolved  in  the  min¬ 
eral  oil  at  a  concentration  C  =  2.3xl0-2  M,  well  above  its  critical  micelle  concen¬ 
tration  (CMC)  ( CMC  =  2.3  xlO-4  M).  At  and  above  the  CMC  the  equilibrium 
tension  7  is  3  mN/m  and  the  equilibrium  surface  concentration  Ycmc  —  3.6  x  10-6 
mole/m2;  dynamic  tension  measurements  indicate  the  desorption  rate  constant, 
kd ,  is  O(10-3  sec-1)  (132).  The  time  scale  for  convection  of  surfactant  along  the 
droplet  a/V,  O(10-1  sec),  is  much  shorter  than  the  desorption  time,  1  /kd,  there- 
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Figure  3.2:  Dielectrophoretic  merging  of  40  /an  radius  droplets  at  the  oil/PDMS 
surface:  Frame  captures  of  the  pairwise  merging  at  time  intervals  of  .12  sec,  flow 
direction  from  bottom  to  top  and  time  (t)  as  function  of  the  measured  edge-to-edge 
scaled  separation  s/a  from  the  images.  The  continuous  line  is  a  fit  for  a  value  for 
the  droplet-wall  drag  coefficient,  a. 


fore  surfactant  collects  at  the  trailing  edge  of  the  droplet,  causing  the  tension  to 

be  larger  at  the  front  than  the  back.  This  Marangoni  gradient  opposes  the  surface 

flow  and  immobilizes  the  interface  because  the  ratio  of  the  characteristic  scale  for 

the  retarding  tension  gradient,  RTTcMc/a ,  to  the  oil  viscous  stress  on  the  droplet 

RTTcmC'  . 


surface,  fi0uV / a ,  the  Marangoni  number  (. Ma  = 


-)  is  O(102).  Hence  the 


pairwise  hydrodynamic  interaction  is  one  of  interfacially  rigid  droplets. 

To  compare  the  data  sets  of  s(t)  to  a  hydrodynamic  model  of  the  merging 
process,  the  applied  held  and  droplet-wall  separation  distance  d  have  to  be  de- 
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termined.  Since  the  PDMS  and  oil  are  dielectric  phases,  they  act  as  capacitors 
in  series,  and  therefore  £  =  V{2 Lpdms(£ou/^pdms)  +  Lc}~ 1  with  £PDMs  and 
eon  the  dielectric  constants  for  the  PDMS  and  oil  (2.65  and  2.18,  respectively), 
and  Lc  and  LPDMS  are,  respectively,  the  length  of  the  channel  and  the  distance 
between  the  electrode  and  the  channel  (6  mm  and  4.5  mm,  respectively).  An  AC 
electric  potential  V  (500  Hz,  5  kV  sine)  is  applied  across  the  electrodes,  result¬ 
ing  in  an  average  held  strength  £  =  365  V/rnm.  While  an  alternating  potential 
is  used  to  prevent  any  residual  charge  accumulation  in  the  PDMS  and  oil,  the 
oscillation  does  not  affect  the  merging  process  since  FPep  is  proportional  to  the 
square  of  the  held,  and  the  oscillation  period  is  much  faster  than  than  time  scale 
for  the  merging.  To  obtain  d,  we  assume  that  at  the  how-focusing  orifice  the 
droplets  detach  symmetrically  from  the  top  and  bottom  walls  of  the  channel  and 
are  therefore  initially  centered  at  the  midplane,  hence  the  initial  separation  d,  is 
h/2  —  a.  Clearly,  each  drop  is  not  precisely  released  at  the  center  of  the  channel, 
but  statistically  above  and  below  the  midplane.  However,  our  data  analysis  will 
average  over  several  data  sets  at  a  nominal  value  of  d  which  should  account  for 
this  statistical  variation.  The  distance  d,  accounting  for  only  the  resistance  of 
the  lower  wall,  is,  for  d/a  «  1,  given  by  (133)  in  {  =  -2ga(p^~poi0r  where 

g  is  the  acceleration  of  gravity.  (Values  for  d/a  are  less  than  0.25,  for  which  this 
expression  is  accurate.  The  equation  for  d  can  be  corrected  to  include  wall-slip 
and  the  effect  of  the  resistance  due  to  the  opposite  channel  wall,  but  these  correc¬ 
tions,  if  included,  can  be  shown  to  be  negligible  in  the  determination  of  d  (134). 

The  droplet  hydrodynamics  is  in  the  Stokes  regime  of  negligible  inertia  (Re  = 

—  =  O(10~4)),  and  the  droplets  remain  spherical  until  the  onset  of  coa- 
lescence  (cf.  Fig.  3.2)  as  the  viscous  forces  are  smaller  than  the  tension  force 
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(capillary  number,  Ca  =  —  =  0(10  2))  and  the  Maxwell  electrical  stresses  are 

7 

smaller  than  the  tension  force  (electric  Bond  number,  Bo  =  -  — - =  0(1CT2)), 

7 

where  eQ  is  the  permittivity  of  free  space.  The  total  fluid  drag  exerted  on  each 

of  the  droplets  as  they  merge  can  be  calculated  from  the  sum  of  the  fluid  drags 

(formulated  as  a  drag  coefficient  /  multiplied  by  67r fj,0ua  and  a  velocity)  in  three 

flow  configurations  (Fig.  3.3):  the  Poiscuille  flow  over  fixed  droplets  (fp),  and 

the  motions,  stationary  in  the  farfield,  of  A  with  B  fixed  or  B  with  A  fixed,  with 

the  later  two  configurations  each  divided  into  a  mutual  approach  (fm)  and  an  in- 

2  Fjjep 

tandem  motion  (/„).  The  total  fluid  drag  balances  Fdep,  Va  —  Vb  = - — . 

6  nnoajm 

The  slip  is  obtained  by  comparison  of  a  theoretical  calculation  of  fm  (a  function 
of  h/a,  s/a,  d/a  and  A/a)  to  an  experimental  value  calculated  through  the  mea¬ 
surement  of  Va  —  Vb  (or  equivalently  s(t))  and  Fdep- 

To  obtain  Fdep,  the  electric  field  in  the  mineral  oil  around  the  merging  (un¬ 
charged)  water  droplets  in  the  microchannel  is  approximated  by  the  bispherical 
harmonic  solution  (135)  for  the  electrostatic  field  in  an  unbounded  insulating  di¬ 
electric  surrounding  a  pair  of  perfectly  conducting  spheres  (zero  net  charge)  due 
to  a  farfield  electric  field  applied  in  the  direction  of  the  line  of  centers  between 
the  spheres,  from  which  Fdep  =  £o£oii0-2E,2<A(s/a).  $s(s/a)  is  an  infinite  series 
function  of  s/a  evaluated  here  with  polynomial  interpolation  (see  Appendix  B). 

The  coefficient  fm  is  calculated  numerically  from  a  COMSOL  (4.2a)  finite 
element  simulation  using  the  experimental  parameters.  The  simulation  is  first 
verified  by  computing  the  drag  coefficient  for  a  single  sphere  in  a  plane-parallel 
channel  moving  parallel  to  the  wall,  denoted  as  a(d/a,  X/a,h/a),  and  compar¬ 
ing  for  A/a  =  0,  to  an  interpolating  formulae  for  a  for  no  slip  obtained  from  a 
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Figure  3.3:  The  drag  coefficient  fm  as  a  function  of  s/a  for  no-slip  (A  =  0)  and 
A  =  1  /rm  for  d  =  1  /jm  (top)  and  d  =  100  nm  (bottom)  for  h  =  100  /xm  and  a  =  40 
/jin. 
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multipole  solution  by  Feuillebois  et  al  (1)  (see  Appendix  B).  Fig.  3.3  shows  fm 
(symbols)  as  a  function  of  s/a  for  no  slip  and  A  =  1  /an,  for  a  =  40  /i m  and  for 
a  separation  distance  d  =  1  /mi  (top)  and  100  nm  (bottom).  As  expected,  the 
closer  the  droplet  pair  to  the  wall,  the  greater  is  the  influence  of  the  slip,  and  it  is 
clear  that  when  d  of  the  order  of  a  few  hundred  nanometers  a  order  ten  percent 
reduction  in  fm  is  achieved  from  the  no-slip  case  for  A  =  1  /mi,  and  this  change  is 
the  basis  of  our  measurement  of  a  order  one  micron  slip  length.  A  slip  coefficient 
can  be  obtained  from  the  experimental  profiles  s(t)  by  comparison  to  the  inte¬ 
gration  of  —  =  Va  —  Vb  =  A  DEP  using  the  numerical  calculations  of  fm  as  a 
dt  07r/i0a/m 

function  of  s.  We  avoid  this  extended  calculation  by  an  approximation  for  fm  as 
the  sum  of  the  drag  on  a  single  (rigid)  droplet  moving  at  a  distance  d  in  a  channel 
of  height  h  with  slip  A,  a(d/a,  A/a,  h/a )  (obtained  by  COMSOL  calculation  and 
independent  of  s),  the  drag  on  a  (rigid)  droplet  pair  mutually  approaching  at  a 
distance  s  from  each  other  in  an  infinite  medium,  (Jeffreys  solution  for  which 
there  is  a  correlation,  fR(^)  =  (l  +  f^}  jl  +  .38e-{fa“+'68}  /63|,  see  Appendix 
B)  and  a  correction  D(s/a)\ 

/Zp'“  =  «(-,-,-)+*{£}+fl{£}  (3.1) 

{a  a  a )  l  a  J  l  a  J 


The  correction  factor  T>(s/a)  accounts  for  the  double  counting  evident  in  f£PProx 
when  s/a  tends  to  infinity  and  should  tend  to  -1;  we  use  D(s/a)  =  —  which 
allows  congruence  of  the  approximate  formulation  over  the  entire  range  of  s/a, 
cf.  Fig.  3.3.  A  theoretical  prediction  for  s(t)  can  be  constructed  by  integrating 
the  force  balance,  using  f%f/prox: 


s/a 


Si/a 


n  f  —  A  |  f  s  ]  |  n  f  g  ) 

X  a1  a 1  a  I _ la  J _ l  a  J 

£0£0il£2S{^}  /S-Kfloii 
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(a) 


(b) 

Figure  3.4:  Hydrodynamic  drag  coefficient  a  as  a  function  of  separation  d/a  for 
droplets  of  radius  37.5  /im  (a)  and  40  /xm  (b)  for  fixed  height  /i=100  //m.  Symbols 
with  error  bars  are  the  experimentally  fitted  coefficients,  the  remaining  symbols 
are  from  numerical  simulation  (the  acc6&panying  dotted  lines  are  a  guide),  and 
the  no-slip  line  is  from  the  Feuillebois  et  al  (1)  correlation. 


3.1  Background 


where  S{  is  the  initial  separation  of  the  droplet  pair.  This  prediction  is  easily  fit 
to  a  data  set  by  adjusting  a  (which  is  independent  of  s),  and  the  fit  is  shown 
for  the  illustrative  data  in  Fig.  3.2.  (In  practice,  to  include  the  forces  on  the 
droplets  clue  to  immediate  neighbors  of  the  train  at  an  assumed  distance  st  from 
the  pair,  we  have  addended  their  dielectrophoretic  dipole  contribution  ( ^  ) 

to  A  and  the  leading  order  droplet-droplet  interaction  1  +  a/(2st)  to  A.)  All 
data  sets  corresponding  to  different  d  and  a  are  fit  in  this  way.  To  correlate 
the  fitted  values  of  a  to  A,  we  first  bin  all  the  data  sets  into  groups  in  which  in 
each  group  the  radii  differ  by  at  most  3  percent.  Each  binned  group  in  diameter 
is  then  further  binned  into  wall  separation  distances  d  which  differ  by  no  more 
than  3  percent.  The  results  are  plotted  as  the  symbols  with  error  bars  (from  the 
standard  deviation  from  the  average)  in  Fig.  3.4  for  a  droplet  radii  bin  a=37.5± 
1.25  prn  and  a  bin  a=40±1.25  pm  as  a  function  of  the  (binned)  values  of  d/a. 
As  each  droplet  radius  bin  corresponds  to  a  fixed  value  of  h/a ,  the  theoretical 
value  of  a  in  these  bins  is  only  a  function  of  d/a  and  A/a.  Plotted  in  Fig.  3.4 
as  symbols  (for  the  two  radii  bins)  are  the  theoretical  values  as  a  function  of  d/a 
for  values  of  A  equal  to  zero  (the  Feuillebois  correlation  (1)),  1  and  4  pm  (from 
COMSOL  calculation).  For  the  PDMS  bottom  channel  wall,  the  comparison  of 
the  theoretical  and  experimental  values  of  a  show  clearly  a  micron-sized  slip,  while 
for  a  bottom  microchannel  wall  made  of  glass,  no-slip  is  obtained  as  expected. 

Our  demonstration  of  0(1  pm)  slip  at  the  interface  of  an  oil  and  a  polymeric 
surface  (PDMS)  which  releases  to  the  surface  air  retained  in  the  material  to  form 
a  lubricating  layer,  may  serve  as  method  for  enabling  giant  slip  without  having  to 
modify  the  surface  with  an  air-sequestering  texture.  These  results  are  particularly 
relevant  to  microfluidics  where  PDMS  is  the  standard  material,  and  the  use  of 
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oil  streams  with  reagent  water  droplets  have  become  a  dominant  lab  on  a  chip 
platform. 
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4 

Electrocoalescence  of 
Water-in-Crude  Oil  Emulsions  in 
Two  Dimensions 

4.1  Background 

Emulsion  stability  is  relevant  to  a  wide  range  of  applications,  including  foods, 
cosmetics,  petroleum  and  other  industrial  processes  (136,  137,  138).  The  funda¬ 
mental  issue  in  the  understanding  of  emulsion  stability  is  the  interaction  between 
droplets  (or  bubbles)  of  the  dispersed  phase,  specifically  the  time  required  for 
them  to  approach  and  coalesce.  This  can  be  understood  in  terms  of  the  time 
required  to  drain  the  continuous  fluid  phase  from  between  adjacent  droplets.  To 
predict  this  time,  knowledge  of  the  forces  between  the  droplets  is  required  to 
make  the  coalescence  time  a  well-defined  hydrodynamics  problem.  Validation  of 
such  a  prediction  is  most  accurately  realized  by  direct  observation  of  individual 
droplet  pair  interactions,  but  this  can  be  difficult  to  achieve  from  experiments 
performed  on  bulk  emulsions. 

Microfluidics  offers  the  ability  to  generate  emulsions  of  monodisperse  droplets 
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and  observe  interactions  between  droplets  on  a  pairwise  basis  (131).  Droplet  sep¬ 
aration  distances  can  be  precisely  measured,  enabling  calculation  of  droplet  inter¬ 
action  forces  on  a  pairwise  basis.  The  origin  of  these  interaction  forces  depends 
on  the  emulsion  system  being  studied  and  its  application.  In  some  instances, 
emulsion  stability  is  desired,  while  other  applications  are  predicated  on  efficient 
separation  of  the  emulsified  phases.  In  these  applications,  an  external  held  is 
often  applied  to  facilitate  droplet  coalescence  and  separation. 

A  primary  example  of  this  field- assisted  separation  is  clectrocoalescence.  The 
importance  of  electrocoalescence  is  primarily  due  to  its  utilization  in  the  oil  in¬ 
dustry,  where  it  is  used  to  separate  water  from  crude  oil  in  a  process  colloquially 
known  as  ’desalting’.  Water  droplets  intentionally  introduced  into  the  crude  oil 
to  extract  salt  are  subsequently  removed  via  application  of  an  electric  field,  which 
polarizes  the  conducting  droplets,  producing  attractive  (and  repulsive)  forces  be¬ 
tween  droplets  that  result  in  coalescence  of  droplet  pairs.  Successive  coalescence 
events  between  droplets  produce  progressively  larger  droplets  that  settle  under 
gravity  into  a  bulk  water  phase  that  can  be  easily  removed.  The  simplest  model  of 
clectrocoalescence  assumes  that  the  water  droplets  behave  as  infinitely  conduct¬ 
ing  spheres  in  a  perfectly  insulating  oil,  and  that  the  interfacial  rheology  of  the 
water-crude  oil  interface  plays  no  role  in  droplet  coalescence.  In  practice,  many 
crude  oils  have  significant  conductivity,  reducing  droplet  polarization  when  the 
electric  field  is  applied  at  0  ffz.  More  importantly,  crude  oils  contain  numerous 
polar  compounds  known  as  asphaltenes,  which  adsorb  at  the  oil-water  interface 
and  act  to  stabilize  the  droplets,  inhibiting  coalescence.  Due  the  unique  nature 
of  the  asphaltene  constituents  in  a  specific  crude  oil,  the  effect  of  these  surface- 
active  compounds  cannot  be  ascertained  independently.  A  valid  experimental 
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investigation  should  therefore  be  conducted  using  the  crude  oil  under  study  to 
accurately  predict  electrocoalescence. 

Numerous  studies  of  various  aspects  of  electrocoalescence  have  been  published. 
Several  examine  the  physics  of  the  coalescence  event  itself  (139,  140,  141).  Oth¬ 
ers  examine  the  behavior  of  droplet  pairs  initially  separated  by  some  distance. 
Bibette  et  al.  measured  the  effect  of  field  strength  and  droplet  separation  dis¬ 
tance  on  the  coalescence  of  microfluidically  generated  surfactant-stabilized  water 
droplet  pairs  in  hexadecane,  mapping  out  a  phase  diagram  which  includes  three 
regimes:  coalescing,  non-coalescing  and  partial  coalescing  (142).  Chiesa  et  al. 
modeled  the  trajectory  of  a  small  droplet  coalescing  into  a  larger  one  as  a  force 
balance  between  the  attractive  electrostatic  and  resistive  hydrodynamic  forces, 
showing  good  agreement  between  experiment  and  theory  (143).  However,  no 
current  study  examines  electrocoalescence  in  two  dimensions,  with  interactions 
between  multiple  droplets  taken  into  account.  Furthermore,  all  of  the  existing 
studies  are  performed  in  model  systems,  which  are  of  limited  applicability  to  wa¬ 
ter  in  crude  oil  emulsions  where  the  heterogeneous  nature  of  the  crude  oil  may 
have  a  dramatic  impact  on  electrocoalescence. 

This  study  reports  results  of  electrocoalescence  experiments  performed  on  two 
dimensional  configurations  of  monodisperse  water  droplets  in  a  crude  oil.  The  ob¬ 
jective  was  to  accurately  predict  electrocoalescence  between  droplet  pairs  based 
on  a  calculation  of  the  pairwise  electrostatic  force  as  determined  by  finite  element 
simulations. 


63 


4.  ELECTROCOALESCENCE  OF  WATER-IN-CRUDE  OIL 
EMULSIONS  IN  TWO  DIMENSIONS 


Figure  4.1:  Schematic  of  PDMS  microchannel  geometry  used  in  experiments. 

4.2  Experimental  Setup 

To  enable  direct  observation  of  water  in  crude  oil  emulsions,  experiments  were 
performed  in  microfluidic  channels  composed  of  PDMS.  The  microchannels  were 
fabricated  using  standard  soft  lithography  techniques  to  produce  two  layers  of 
PDMS,  one  containing  the  fluidic  channel  and  a  second,  flat  layer  to  seal  the 
channel.  The  two  layers  were  bonded  together  following  exposure  to  an  oxygen 
plasma  and  mounted  on  a  standard  glass  microscope  slide.  Access  ports  cored 
into  the  top  layer  using  a  biopsy  punch  allowed  introduction  fluids  of  interest  via 
polyethylene  tubing  connected  to  syringe  pumps.  To  generate  electric  fields  in 
the  microchannel,  planar  electrodes  of  aluminum  were  inserted  into  the  PDMS 
perpendicular  to  the  lateral  plane  of  the  microchannel  and  sited  externally  from 
the  fluid  (Figure  4.1).  The  electrodes  were  connected  to  an  amplifier  (Trek)  con¬ 
trolled  by  a  frequency  generator  (Agilent). 

The  first  experiment  demonstrated  the  ability  to  directly  observe  electrocoa¬ 
lescence  of  a  water  in  crude  oil  emulsion  in  a  microfluidic  channel.  The  emulsion 
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Figure  4.2:  Time  sequence  of  electrocoalescence  in  water-in-crude  oil  emulsion. 

was  generated  by  mixing  water  and  crude  oil  in  a  prescribed  volumetric  ratio  in  a 
blender.  The  emulsion  was  then  loaded  into  a  glass  syringe  and  connected  to  the 
microfluidic  device,  which  was  mounted  on  an  inverted  brightfield  microscope. 
The  emulsion  was  flowed  through  the  microchannel  and  a  uniform  electric  field 
E  was  applied  across  the  channel.  Images  of  the  emulsion  were  captured  using  a 
high-speed  camera.  Individual  coalescence  events  between  water  droplets  in  crude 
oil  were  observed  (Figure  4.2),  as  was  the  change  in  droplet  size  distribution  of 
the  emulsion  over  time,  demonstrating  successful  observation  of  the  electrocoa¬ 
lescence  in  crude  oil. 

Due  to  the  polydisperse  nature  of  a  conventional  emulsion,  direct  observation 
of  the  individual  electrocoalescence  events  between  droplets  is  of  limited  utility. 
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Ideally,  emulsion  droplets  should  be  monodisperse  and  have  a  known  orienta¬ 
tion  to  one  another,  enabling  quantification  of  the  forces  between  them.  This  is 
achieved  by  utilizing  a  microfluidic  technique  known  as  flow  focusing.  By  intro¬ 
ducing  the  water  and  crude  oil  phases  separately  and  flowing  them  through  an 
orifice,  monodisperse  droplets  of  water  in  crude  can  be  formed.  This  allows  the 
subsequent  manipulation  of  the  droplets  by  the  electric  field  to  be  accurately  mod¬ 
eled  because  both  the  droplet  size  and  the  separation  distances  between  droplets 
are  known.  To  generate  two-dimensional  configurations  of  droplets,  the  one  di¬ 
mensional  droplet  train  formed  immediately  downstream  from  the  flow-focusing 
orifice  is  directed  into  a  large  channel  where  the  droplet  velocities  decrease  and 
the  droplets  are  arranged  in  an  arbitrary  configuration  (Figure  4.1).  The  electric 
field  is  applied  across  the  channel  and  the  droplet  coalescence  events  are  recorded 
using  a  high-speed  camera.  To  prevent  the  electrical  conductivity  of  the  crude 
oil  from  reducing  the  electric  field  over  time,  the  electric  field  is  applied  at  a 
frequency  of  500  Hz. 


4.3  Droplet  Force  Calculation 


To  predict  electrocoalescence  between  the  water  droplets  in  the  two-dimensional 

configuration  upon  application  of  the  electric  field,  the  electrostatic  forces  between 

the  droplet  pairs  must  be  calculated.  The  magnitude  of  the  applied  electric  field 

is  less  than  that  required  to  distort  the  shapes  of  the  droplets  (the  electric  Bond 
e0eOuE2a  . 


number,  Bo  = 


7 


-,  is  less  than  one),  therefore  the  droplets  are  modeled 


as  monodisperse  spheres  with  known  diameter  and  position.  The  period  of  the 
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Figure  4.3:  Two  dimensional  configuration  of  water  droplets  in  crude  oil  with 
COMSOL  model. 

applied  field  is  much  greater  than  the  polarization  time  of  the  droplets,  r  = 
£h2o/°'H20-  This  produces  an  electric  field  distribution  around  the  droplets  that 
is  quasi-static  and  can  be  calculated  from  the  solution  to  the  Laplace  equation 
subject  to  the  boundary  conditions  at  the  droplet  surfaces: 

V2V  =  0  (4.1) 

Von  =  VH2o  (4.2) 

n-(&OiiEoii  —  CF/2oE_ff2o)  =  0  (4.3) 

Although  the  external  field  is  uniform,  the  field  around  the  droplets  is  non-uniform 
due  to  their  polarization,  with  larger  held  strengths  at  the  two  poles  aligned  with 
the  external  held  and  lower  held  strengths  at  the  two  poles  aligned  perpendicular 
to  the  external  held.  In  the  case  of  a  single  droplet,  the  magnitudes  of  the  held 
strengths  at  these  poles  will  be  equal,  resulting  in  no  net  force  on  the  droplet  and 
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therefore  no  motion.  However,  droplets  in  proximity  to  one  another  will  experi¬ 
ence  non-symmetric  local  fields  and  will  consequently  experience  net  attractive 
or  repulsive  forces  depending  on  their  orientation  relative  to  the  external  held. 
In  the  simple  case  of  a  single  droplet  pair,  this  net  force  has  been  theoretically 
calculated  using  bispherical  coordinates  (135).  In  the  case  of  a  configuration  of 
multiple  droplets,  finite  element  simulations  are  used.  These  simulations  solve 
the  Laplace  equation  for  each  configuration  of  droplets  to  calculate  the  electric 
held  E  around  each  droplet  surface.  This  held  is  used  to  calculate  the  net  force 
on  each  droplet. 

F Dr0piet  =  <J)n-TdA  -  j)p-ndA  (4.4) 

where  n  is  the  surface  normal  and  T  is  the  Maxwell  stress  tensor: 


T  =  eEE  -  -Ea« 


(4.5) 


Applying  the  Maxwell  equations  V-eE  =  0  and  VxE  =  0  and  a  vector  identity 
simplifies  the  divergence  of  T.  which  is  equal  to  the  pressure  gradient: 


_  £ _  e°V7lTT'|2 


V|E|  =  Vp 


(4.6) 


E  has  no  tangential  component,  therefore  |E | 2  =  En 2  and  the  net  force  on  the 
droplet  can  be  calculated: 


F  Droplet  s/^  ~  En  (LA 


(4.7) 


This  net  force  is  expected  to  predict  the  droplet  motion  upon  application  of  the 
electric  held.  However,  the  close  proximity  of  the  droplets  prior  to  the  application 
of  the  held  means  that  no  significant  motion  between  droplets  is  possible  prior  to 
coalescence.  Therefore,  the  net  force  on  each  droplet  must  be  deconstructed  to 


4.3  Droplet  Force  Calculation 


identify  the  component  of  the  force  due  to  each  adjacent  droplet.  In  other  words, 
the  forces  on  the  droplet  must  be  calculated  on  a  pairwise  basis  to  identify  which 
droplet  pairs  experience  net  attractive  forces  and  are  therefore  likely  to  coalesce. 
This  is  done  by  simple  vector  projection  using  the  net  force  of  each  droplet  and 
the  center  to  center  vector  between  each  droplet  pair: 


F  = 

±  pair 


|F2  —  Fi | •  |r2  -  rij 

r2  -  ri 


(4.8) 


Because  each  droplet  only  has  a  certain  number  of  nearest  neighbors  which  block 
its  potential  coalescence  with  other  droplets,  an  algorithm  is  employed  to  identify 
which  droplet  pairs  are  not  obstructed  by  third  droplets  and  are  therefore  poten¬ 
tially  able  to  coalescence.  Using  these  calculations,  the  resulting  pairwise  droplet 
forces  are  plotted  as  a  function  of  the  separation  distance  between  the  droplets 
in  the  pair. 

The  coalescence  of  individual  droplet  pairs  proceeds  upon  application  of  the 
external  electric  held.  Due  to  the  fact  that  the  droplet  configuration  changes 
subsequent  to  initial  coalescence  events,  only  droplets  pairs  coalescing  within  the 
first  0.1  seconds,  which  are  therefore  assumed  to  be  coalescing  due  to  the  forces 
calculated  from  the  initial  droplet  configuration,  are  considered.  These  coalesc¬ 
ing  pairs  are  identified  on  the  plot  of  pairwise  droplet  forces  versus  separation 
distance  to  determine  if  the  calculated  forces  accurately  predict  coalescence.  The 
results  are  shown  in  Figure  4.3  for  a  representative  data  set.  As  expected,  droplet 
pairs  with  large  negative  (attractive)  forces  between  them  coalesce,  while  droplet 
pairs  with  positive  (repulsive)  forces  between  them  do  not  coalesce.  At  larger 
separation  distances,  droplet  pairs  do  not  experience  large  forces  and  therefore 
do  not  coalesce  due  to  the  initial  droplet  configuration.  As  shown  in  Figure  4.4, 
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Figure  4.4:  Plot  of  electrostatic  forces  between  droplet  pairs  versus  normalized 
separation  distance. 
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several  droplet  pairs  with  negative  forces  do  not  coalesce.  This  can  be  explained 
by  the  fact  that  one  droplet  in  the  pair  is  also  half  of  a  pair  with  a  larger  nega¬ 
tive  force,  and  therefore  coalesces  with  the  droplet  with  which  it  has  the  larger 
attractive  force,  as  expected. 

In  conclusion,  we  have  demonstrated  accurate  prediction  of  pairwise  electro¬ 
coalescence  of  water  droplets  emulsified  in  a  crude  oil.  The  unique  ability  of 
microfluidic  geometries  to  generate  monodisperse  emulsion  droplets  and  observe 
them  in  a  macroscopically  opaque  continuous  phase  such  as  crude  oil  should  en¬ 
courage  further  work  using  microfluidics  to  study  emulsion  stability. 
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5 

Future  Work 


This  proposal  outlines  a  new  method  of  measuring  the  rheological  properties  of 
liquid- liquid  interfaces  (e.g.  oil/water  interfaces)  with  adsorbed  surfactant  mono- 
layers  using  a  microfluidic  geometry  via  application  of  an  electric  field.  Monodis- 
perse  water  droplets  formed  in  oil  in  a  microchannel  with  an  adsorbed  monolayer 
at  their  interface  are  subjected  to  uniform,  oscillating  electric  fields  that  dis¬ 
tort  their  shape  according  to  the  balance  between  surface  stresses  and  electrical 
stresses.  By  evaluating  the  change  in  the  shape  of  the  droplets  with  time,  the 
surface  tension  and  the  surface  rheological  dilatational  viscosity  due  to  the  sur¬ 
face  monolayer  can  be  measured.  This  technique  can  be  used  to  study  the  surface 
rheology  of  surface  active  species  at  liquid-liquid  interfaces  in  many  applications, 
including  protein  therapeutics  to  understand  the  relationship  between  adsorp¬ 
tion  and  aggregation.  It  offers  significant  advantages  over  existing  measurement 
techniques  due  to  its  versatility  in  measuring  interfaces  between  two  phases  of 
comparable  densities,  interfaces  with  an  opaque  continuous  phase  and  interfaces 
with  extremely  low  surface  tension  values.  The  high-throughput  characteristics 
of  microfluidic  systems  allow  the  rapid  generation  of  high  number  statistics  and 
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eliminate  the  potential  of  contamination,  while  the  small  length  scales  provide 
the  opportunity  to  study  sensitive  measures  of  the  surface  rheology. 

5.1  Background 

Our  interest  in  the  surface  rheology  of  surfactant-laden  oil-water  interfaces  derives 
from  the  importance  of  this  rheology  in  emulsions,  and  also  in  evaluating  protein 
interactions  (as  will  be  discussed  in  Proposed  Research).  Emulsions  are  disper¬ 
sions  of  two  immiscible  liquids,  typically  water  and  oil,  in  which  droplets  of  one 
liquid  phase  are  dispersed  in  a  continuous  phase  of  the  second  liquid.  Emulsions 
are  present  in  many  applications  where  they  are  either  introduced  intentionally 
(pharmaceuticals,  foods,  cosmetics)  or  where  they  occur  naturally  (petroleum) 
(136,  137,  138).  The  stability  of  the  emulsion  to  coalescence  of  the  droplets,  and 
separation  of  the  phases,  is  the  central  issue  in  all  applications. 

While  a  stable  emulsion  is  desirable  in  some  applications  (drug  delivery,  food 
processing),  other  applications  focus  on  breaking  the  emulsion  to  achieve  a  phase 
separation.  The  stability  of  an  emulsion  can  be  controlled  by  surface  active 
components,  surfactants,  which  are  amphiphilic  molecules  with  polar  and  non¬ 
polar  groups.  Surfactants  adsorb  at  the  liquid-liquid  (e.g.  oil/water)  interface 
to  form  a  monolayer,  straddling  the  surface  with  the  polar  group  in  the  water 
phase  and  the  non-polar  group  in  the  oil  phase  (see  for  example  the  monographs 
(144,  145,  146,  147,  148)).  Surfactant  adsorption  and  monolayer  formation  change 
the  properties  of  the  fluid  interface  and  thereby  affect  the  stability  of  the  emul¬ 
sion,  particularly  its  lifetime  (i.e.  how  long  the  droplets  of  the  emulsion  remain 
dispersed),  with  dramatic  impact  on  processing  applications  (149,  150,  151). 
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The  static  interface  between  the  two  fluids  is  characterized  by  the  equilibrium 
interfacial  tension,  a.  When  surfactants  are  present  and  adsorb  at  the  interface, 
adsorption  lowers  the  equilibrium  tension  from  the  value  of  the  clean  interface,  ac, 
to  a  value  that  is  a  function  of  the  surface  concentration  T  of  the  amphiphiles  at 
the  interface,  i.e.  cr(T),  the  surface  equation  of  state.  When  the  phases  bounding 
the  interface  are  in  motion,  the  interface  can  stretch  (dilate),  or  a  shear  flow  can 
be  set  up  in  the  plane  of  the  interface.  At  surfactant-laden  liquid-liquid  interfaces 
in  motion,  the  interfacial  shear  and  dilation  create  stresses  in  the  plane  of  the 
surface  apart  from  either  the  interfacial  tension  (152,  153)  or  gradients  in  the 
interfacial  tension  due  to  gradients  in  the  surface  concentration  of  surfactant  set 
up  by  the  surface  flow  (Marangoni  stresses).  In  the  simplest  description  (a  New¬ 
tonian  surface  fluid),  these  stresses  can  be  related  to  the  surface  flow  by  surface 
shear  (/xs)  and  dilatational  (ks)  viscosity  coefficients,  which  are  functions  of  the 
surface  concentration.  The  interfacial  stress  at  a  moving  interface  is  therefore 
characterized  by  the  equation  of  state,  which  dictates  the  Marangoni  gradients, 
and  the  viscosity  coefficients,  which  govern  the  surface  rheology. 

Many  techniques  exist  to  measure  the  equilibrium  interfacial  tension  as  a  func¬ 
tion  of  the  bulk  concentration  of  surfactant,  and  thereby  the  surfactant  equation 
of  state  (e.g.  capillary  rise,  Wilhclmy  plate,  static  drop/bubble  shape  analy¬ 
sis,  see  for  example,  (154,  155,  156)).  Measurement  of  the  interfacial  viscosities 
is  more  difficult,  particularly  the  dilatational  viscosity,  which  involves  setting  a 
dilating  flow  along  a  surface  and  measuring  the  interfacial  stress  (157).  The  stan¬ 
dard  technique  for  measuring  the  dilatational  viscosity  is  the  oscillating  pendant 
drop  technique  (158,  159).  This  technique  uses  a  syringe  to  produce  a  droplet 
(or  bubble)  of  one  phase  in  another  and  analyzes  a  2D  silhouette  of  the  axisym- 
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metric  shape  as  a  force  balance  between  buoyancy  and  surface  tension  to  obtain 
the  surface  tension.  The  volume  of  the  pendant  drop  can  be  perturbed  by  small 
amounts  at  low  frequencies  to  set  up  a  dilating  flow.  By  simultaneously  measur¬ 
ing  the  tension  and  the  area  expansion  rate,  the  surface  dilatational  viscosity  can 
be  obtained. 

While  the  pendant  drop  technique  is  extremely  useful,  it  has  inherent  dis¬ 
advantages.  The  precision  of  the  measurement  is  susceptible  to  contamination, 
particularly  when  dynamic  measurements  are  being  made  over  long  periods  of 
time.  As  the  method  relies  on  a  buoyancy  force  between  the  droplet  and  bulk 
phase,  it  cannot  be  used  when  the  two  phases  have  comparable  densities,  a  serious 
limitation  in  some  systems.  Additionally,  extremely  low  surface  tension  values 
(<  5  mN/m)  affect  the  stability  of  the  droplet,  making  accurate  measurements 
difficult.  Perhaps  most  importantly,  however,  the  pendant  drop  technique  cannot 
be  used  when  both  phases  are  opaque  to  visible  light. 

5.2  Proposed  Research 

We  propose  to  develop  a  microfluidic  tensiomenter  to  measure  the  dilatational 
viscosity  of  a  liquid-liquid  interface.  Microfluidic  systems  can  be  used  to  generate 
and  study  emulsion  droplets.  The  flow  focusing  technique  allows  the  creation  of 
highly  monodisperse,  surfactant-stabilized  droplets  and  has  been  demonstrated  in 
a  number  of  applications  (131,  160).  These  include  several  techniques  for  measur¬ 
ing  surface  rheology  (161,  162).  We  propose  to  use  a  microfluidic  flow-focusing 
geometry  to  develop  an  interfacial  rheometer  for  measuring  the  surface  dilata¬ 
tional  viscosity  by  expanding  and  contracting  water  droplets  in  a  dielectric  oil 
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using  an  electric  field  (Figure  5.1).  Trains  of  highly  monodisperse  water  droplets 
(d  ~  50-100  /mi)  are  dispersed  in  a  continuous  oil  phase  as  the  two  phases  flow 
through  an  orifice  under  pressure-driven  flow.  The  droplet  train  enters  a  chan¬ 
nel  located  between  two  externally  sited,  parallel  electrodes.  The  electrodes  are 
connected  to  a  frequency  generator  and  amplifier,  and  generate  high  strength  (« 
5  kV/mm)  electric  fields  to  elongate  the  shape  of  the  droplets,  which  are  imaged 
using  a  high-speed  camera.  To  eliminate  the  effect  of  the  shear  flow  on  the  droplet 
shape,  the  flow  can  be  stopped  via  pressure  equilibration  of  the  microchannel  in¬ 
let  and  exit  or  via  the  incorporation  of  a  microfluidic  valve  within  the  channel 
itself. 


Aqueous 

Phase 


microfluidic  cell 


Figure  5.1:  Schematic  of  interfacial  rheometer  using  oscillating  electric  field. 


The  static  shape  of  a  water  droplet  in  an  oil  (assumed  to  be  a  dielectric) 
elongated  by  a  uniform  electric  field  is  determined  by  the  balance  between  the 
capillary  pressure  Pc  =  a / a  and  the  electric  Maxwell  stress  Em  =  eo €quE2  ex- 
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erted  at  the  interface  where  eo  is  the  electrical  vacuum  permittivity,  soil  is  the 

relative  permittivity  of  the  oil,  E  is  the  magnitude  of  the  applied  field  and  a  is 

the  spherical  radius  of  the  droplet  before  deformation.  The  steady  elongation  of 

conducting  and  dielectric  drops  in  a  dielectric  liquid  has  been  studied  extensively 

in  the  literature,  beginning  with  the  analysis  of  Taylor  (163,  164,  165,  166,  167, 

168,  169,  170,  171,  172).  This  elongation  is  analogous  to  a  water  droplet  at  a 

needle  tip  in  oil  which  is  hanging  and  elongated  by  gravity  and  whose  shape  is 

determined  by  a  balance  of  capillary  pressure  and  gravity  as  given  by  the  Young- 

Laplace  equation.  For  the  case  of  the  electric  distortion,  this  balance  is  scaled 

by  the  electric  Bond  number,  Bo  =  —  — - -.  When  the  surface  tension  is  high, 

cr 

the  droplet  will  maintain  its  spherical  shape.  When  the  surface  tension  is  low, 
the  electric  field  will  distort  its  shape  in  the  direction  of  the  field. 

By  applying  a  prescribed  DC  electric  field,  the  equilibrium  shape  of  the  droplet 
can  be  used  to  calculate  the  surface  tension.  The  equation  for  the  axisymmetric 
shape  of  the  surface,  which  is  a  modified  Young-Laplace  equation  with  gravity 
replaced  by  the  electrical  stress,  can  be  solved  to  obtain  the  shape  as  a  function 
of  Bo.  This  calculation  is  more  difficult  than  the  shape  calculation  in  the  case 
of  gravity  because  the  electric  field  on  the  surface  must  be  determined  in  order 
to  calculate  the  Maxwell  stresses.  This  held  is  computed  numerically  by  solving 
the  Laplace  equation  for  the  electric  held  in  the  dielectric  space  around  the  drop 
subject  to  the  electric  held  conditions  on  the  drop  surface  and  the  imposition  of 
a  uniform  held  at  infinity  as  discussed  in  Chapter  4.  As  this  space  is  defined  by 
the  contour  of  the  drop  boundary,  the  calculations  of  the  shape  and  the  held  are 
coupled  and  must  be  solved  numerically. 

We  propose  to  use  COMSOL  finite  element  simulations  to  undertake  these 
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calculations,  and  preliminary  results  have  been  obtained  (see  below).  By  mea¬ 
surement  of  the  shape  through  optical  microscopy  of  the  transparent  microfluidic 
cell,  and  by  comparison  to  the  numerical  solution,  the  electric  Bond  number  Bo 
can  be  obtained.  Since  the  magnitude  of  the  applied  electric  field  E  is  known,  as 
is  the  radius  a,  the  surface  tension  cr  of  the  interface  can  be  calculated. 

To  interrogate  the  interfacial  rheology  of  the  droplet,  the  electric  field  can 
be  oscillated  at  low  frequencies  in  a  manner  directly  analogous  to  the  pressure- 
driven  oscillations  of  the  pendant  drop  technique.  The  elongation  of  the  droplet 
stretches  the  interface,  and  is  resisted  by  the  dilatational  viscosity,  which  adds  to 
the  tension  of  the  drop.  The  amplitude  of  oscillation  must  be  sufficiently  large 
to  produce  an  observable  change  in  droplet  shape,  but  small  enough  so  that  the 
convective  effects  induced  by  the  change  in  shape,  which  are  proportional  to  the 
droplet  size  and  the  oscillation  frequency,  can  be  neglected.  To  determine  the 
optimum  electric  field  amplitude  and  frequency  for  a  given  droplet  size,  mea¬ 
surements  can  be  performed  on  a  clean  droplet  interface  with  a  known  surface 
tension  to  confirm  that  convective  effects  are  not  biasing  the  measurement  (this 
same  procedure  is  employed  with  the  pendant  drop  apparatus). 

An  additional  consideration  is  the  possibility  of  producing  of  Marangoni  forces 
due  to  surface  concentration  gradients  induced  by  the  change  in  droplet  shape 
or  the  migration  of  charged  surfactant  due  to  the  electric  field.  As  in  the  oscil¬ 
lating  pendant  drop,  the  surface  concentration  is  assumed  to  be  sufficiently  large 
to  minimize  the  Marangoni  effects,  and  the  fact  that  the  continuous  oil  phase  is 
assumed  to  be  a  perfect  dielectric  means  that  charged  surfactants  can  only  be 
present  in  the  droplet  water  phase,  where  the  electric  field  is  zero  and  therefore 
charge  migration  is  not  possible. 


79 


5.  FUTURE  WORK 


Under  these  conditions,  the  static  drop  shape  equations  are  valid  and  the  net 
isotropic  tension  (7)  of  the  droplet  can  be  computed  from  each  of  the  oscillating 
shapes  to  obtain  the  oscillating  tension  as  a  function  of  time.  From  the  Newto¬ 
nian  surface  equation  of  state,  the  dilatational  viscosity  ks  can  be  recomputed 
from  the  net  tension  7  and  the  measured  area  expansion  by  neglecting  shear 
contributions  for  the  net  isotropic  tension  7: 

1  dA 

7  =  o-(r)  +  ksVs-v  =  cr(r)  +  — —  (5.1) 

where  v  is  the  velocity  of  the  surface. 

A  microfluidic  system  has  several  advantages  over  the  pendant  drop  tech¬ 
nique.  The  droplet  interfaces  are  continuously  generated  via  flow-focusing,  and 
are  therefore  not  subject  to  contamination.  The  force  balance  does  not  require  a 
significant  buoyancy  force,  allowing  for  investigation  of  liquids  with  comparable 
densities.  The  laminar  flow  regime  in  which  the  droplets  are  formed  allows  gener¬ 
ation  of  droplet  interfaces  with  extremely  low  surface  tensions.  By  controlling  the 
thickness  of  the  microfluidic  channels  in  which  the  experiment  is  performed,  liq¬ 
uids  that  would  appear  opaque  in  larger  path  lengths  are  semi-translucent  (even 
crude  oils).  Finally,  the  advantages  inherent  to  all  microfluidic  systems,  conserved 
analyte  and  high-throughput,  allow  for  the  enhanced  precision  of  high-n  statistics 
without  increases  in  cost  or  time. 

To  test  the  model,  we  will  study  the  surface  rheology  of  emulsions  of  pro¬ 
tein  solutions.  Proteins  adsorbing  at  oil-water  interfaces  have  been  shown  to 
demonstrate  dilatational  properties,  which  may  correlate  to  aggregation  behavior 
(173,  174,  175).  For  this  model,  we  will  use  alpha- Chymotripsinogen  (aCgn),  a 
commercially  available  and  inexpensive  alpha-beta  protein  that  is  known  to  dena- 
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ture  and  aggregate  upon  exposure  to  a  hydrophobic  surface.  Understanding  such 
aggregation  behavior  is  important  for  the  growing  held  of  protein  therapeutics, 
and  may  also  enhance  understanding  of  aggregation- mediated  diseases  (176,  177). 

5.3  Preliminary  Results 

A  key  step  in  the  calculation  of  the  surface  viscosity  from  the  shapes  of  the  water 
droplets  in  oil  driven  by  the  oscillating  electric  held  is  the  computation  of  the 
shape  from  a  modified  Young-Laplace  equation  (with  gravity  replaced  by  the 
electric  held  stresses).  This  calculation  is  coupled  to  the  solution  of  the  Laplace 
equation  for  the  held  in  the  oil  outside  the  drop  subject  to  the  held  conditions  on 
the  drop  surface.  Figure  5.2  shows  a  simple  example  of  the  steady  and  unsteady 
shapes  of  droplets  subject  to  a  uniform  electric  held  from  a  COMSOL  calculation 
in  which  the  interface  is  described  through  a  level  set  formulation,  and  the  Laplace 
equation  for  the  electric  held  is  solved  by  a  finite  element  algorithm  coupled  to  the 
level  set  formulation.  The  electric  Bond  number  Bo  is  larger  than  one,  causing 
the  droplet  to  deform  from  its  spherical  shape  to  a  new  equilibrium  shape.  As 
the  surface  tension  decreases  from  30  mN/m  to  15  mN/m,  Bo  increases  and  the 
deformation  becomes  larger. 
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Figure  5.2:  COMSOL  calculation  of  droplet  in  oil  subject  to  uniform  electric 
field  ( eon  =  2.5,  E  =  500  V/mm,  a  =  1.5  mm),  (a)  a  =  30  mN/m.  (b)  a  =  15 
mN/m.  Spherical  droplets  deform  due  to  the  applied  held  to  reach  a  non-spherical 
equilibrium  shape.  Note  the  longer  time  required  for  the  lower  tension  interface  to 
reach  equilibrium. 
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The  following  data  represent  the  complete  set  of  normalized  surface  concentra¬ 
tion  binding  curves  for  all  parameter  values  used  in  the  finite  element  simulations 
performed  in  Chapter  2.  Kinetically  limited  binding  curves  are  shown  for  com¬ 
parison. 


Varied  Pe,  Da  =  1,  e  =  0.016 


Figure  A.l:  Normalized  binding  curves  showing  effect  of  Pe  at  Da  =  1 

The  following  data  represent  the  complete  set  of  normalized  concentration 
profiles  for  all  parameter  values  used  in  the  finite  element  simulations  performed 
in  Chapter  2. 
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Varied  Pe,  Da  =  10,  e  =  0.016 


Figure  A. 2:  Normalized  binding  curves  showing  effect  of  Pe  at  Da  =  10 


Varied  Pe,  Da  =  100,  e  =  0.016 


Figure  A. 3:  Normalized  binding  curves  showing  effect  of  Pe  at  Da  =  100 


84 


Varied  Da,  Pe  =  10000,  e  =  0.016 


Figure  A. 4:  Normalized  binding  curves  showing  effect  of  Da  at  Pe  =  10000 


Varied  Da,  Pe  =  1000,  e  =  0.016 


Figure  A. 5:  Normalized  binding  curves  showing  effect  of  Da  at  Pe  =  1000 
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Varied  Da,  Pe  =  10,  e  =  0.016 


Figure  A. 6:  Normalized  binding  curves  showing  effect  of  Da  at  Pe  =  10 


Figure  A. 7:  Cross-section  of  concentration  profile  in  microchannel  for  patch  sur¬ 
face  at  x  =  0,  Pe  =  100,  Da  =  1 


Figure  A. 8:  Cross-section  of  concentration  profile  in  microchannel  for  bead  sur¬ 
face  at  x  =  0,  Pe  =  100,  Da  =  1 


Figure  A. 9:  Cross-section  of  concentration  profile  in  microchannel  for  patch  sur¬ 
face  at  x  =  0,  Pe  =  100,  Da  =  10 


Figure  A. 10:  Cross-section  of  concentration  profile  in  microchannel  for  bead 
surface  at  x  =  0,  Pe  =  100,  Da  =  10 
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Figure  A. 11:  Cross-section  of  concentration  profile  in  microchannel  for  patch 
surface  at  x  =  0,  Pe  =  100,  Da  =  100 


Figure  A. 12:  Cross-section  of  concentration  profile  in  microchannel  for  bead 
surface  at  x  =  0,  Pe  =  100,  Da  =  100 


Figure  A. 13:  Cross-section  of  concentration  profile  in  microchannel  for  patch 
surface  at  x  =  0,  Pe  =  1000,  Da  =  1 


Figure  A. 14:  Cross-section  of  concentration  profile  in  microchannel  for  bead 
surface  at  x  =  0,  Pe  =  1000,  Da  =  1 


Figure  A.  15:  Cross-section  of  concentration  profile  in  microchannel  for  patch 
surface  at  x  =  0,  Pe  =  1000,  Da  =  10 


Figure  A.  16:  Cross-section  of  concentration  profile  in  microchannel  for  patch 
surface  at  x  =  0,  Pe  =  1000,  Da  =  10 
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Figure  A.  17:  Cross-section  of  concentration  profile  in  microchannel  for  bead 
surface  at  x  =  0,  Pe  =  1000,  Da  =  10 


Figure  A. 18:  Cross-section  of  concentration  profile  in  microchannel  for  patch 
surface  at  x  =  0,  Pe  =  1000,  Da  =  100 
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Figure  A.  19:  Time  sequence  of  fluorescent  micrographs  measuring  binding  of 
NeutrAvidin- Texas  Red  (C  =  4.16  E  gfj)  M )  to  biotin- functionalized  glass  mi¬ 
crobeads  (T  =  5.5  E  -9  M)  at  Pe  =  5600. 
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B.l  Materials 

B.1.1  Aqueous  Phase 

DI  Water:  From  Millipore  ultrafiltration  unit,  conductivity,  18  Mil  cnr1,  and 
from  handbook  values  pw=  103  kg  m-3,  and  fiw—  10-3  kg  m_1s_1  at  20  °C. 
Mineral  Oil:  Sigma-Aldrich,  eoii  =  2.18  (parallel  plate  capacitor  measurement), 
PoU  =  8.75  x  102  kg  m-3  (manufacturer’s  data  sheet)  at  20  °C.  The  oil  viscosity 
Pen  is  temperature  sensitive  in  the  range  in  which  the  experiments  were  under¬ 
taken  (20-23°C),  and  was  measured  (with  the  surfactant,  Span  80,  C=2.3xl0-2 
M)  by  a  temperature-controlled  rheometer  (see  data  below,  Figure  B.l)  and  fit  to 
a  correlation  p{T)  =  2.20  x  10~5T3  -  1.63  x  10"3T2  +  3.58  x  10"2T  -  0.139.  A 
thermocouple  in  the  PDMS  measured  the  temperature  for  each  experiment,  and 
in  the  fitting  for  the  drag  coefficient  a  the  viscosity  at  the  temperature  measured 
for  the  experiment  was  used  by  calculation  from  the  correlation. 
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Figure  B.l:  Viscosity  of  mineral  oil  as  a  function  of  temperature. 


B.1.2  Microfluidic  Cell  Fabrication 

The  microfluidic  cell  is  fabricated  via  the  methods  of  soft  lithography  (178)  from 
two  slabs  of  polymerized  and  cured  PDMS,  one  slab  containing  the  inscribed 
channels  and  a  second,  flat  slab  bonded  to  the  first  to  enclose  the  channels.  The 
PDMS  prepolymer  and  curing  agent  are  obtained  as  a  kit,  Sylgard  184  (Dow 
Corning).  The  top  half  of  the  cell  containing  the  inscribed  fluidic  circuitry  is 
molded  from  a  master  fabricated  by  photolithography  using  a  transparency  mask 
and  a  negative  tone  SU-8  photoresist  and  developer  (Microchem),  and  the  lower 
half  is  a  cast  layer  of  PDMS  which  is  bonded  to  the  top  half  by  plasma  treating 
the  surfaces  and  then  conformally  contacting  the  halves.  Dielectric  constant  of 
cured  elastomer  is  2.65  (manufacturer’s  material  data  sheet).  Fluids  are  delivered 
to  the  cell  from  syringe  pumps  (Harvard  PHD)  via  polyethylene  tubing  (1.5  mm 
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ID)  through  punch-hole  entry  ports. 

B.2  Force  Expressions 

B.2.1  Interpolation  Formula  for  the  Electrostatic  Force 
Between  Conducting  Spheres  as  a  Function  of  Sphere- 
Sphere  Separation  Distance 

The  calculation  of  Davis  (135)  for  the  dielectric  force  on  a  pair  of  zero  charge 
perfectly  conducting  spheres  in  an  unbounded  insulating  dielectric  subject  to  a 
uniform  electric  field  £  at  infinity  applied  along  the  line  of  centers  of  the  spheres 
is  given  as  an  infinite  series  solution  of  the  edge-to-edge  separation  distance  s  di¬ 
vided  by  a,  from  which  the  dielectrophoretic  force  between  the  spheres  is  obtained 
as  Fdep  =  £0£0ii8,2a2<^s(s/a)  where  9( s/a )  is  an  infinite  series.  We  approximate 
the  series  by  the  polynomial  Q/s/a)  =  —2.02  x  10^' (s/a)3  +  7.33  x  10_7(s/a)2  —  9.99  x  10_7(s( 
and  a  comparison  of  this  interpolation  formulae  and  the  Davis  expression  is  given 
in  Figure  B.2. 

B.2. 2  Analytical  Solution  for  the  Hydrodynamic  Drag  Force 
Due  to  Approaching  Spheres  as  a  Function  of  Sphere- 
Sphere  Separation  Distance  in  an  Infinite  Medium, 

fR(s/a) 

The  solution  for  the  hydrodynamic  drag  exerted  on  a  pair  of  two  mutually  ap¬ 
proaching  (interfacially  rigid  or  solid)  spheres  in  an  infinite  medium  is  the  Stimson 
and  Jeffrey  solution  in  bispherical  coordinates  and  is  expressed  as  an  infinite  se¬ 
ries  and  the  drag  coefficient  fR(s/a)  is  expressed  as  an  infinite  series  (see  eq.  B.l); 
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Figure  B.2:  Normalized  electrostatic  force  as  a  function  of  normalized  separa¬ 
tion  distance  as  given  by  the  exact  bispherical  calculation  and  the  interpolating 
equation. 


for  numerical  calculation  we  use  the  interpolation  as  given  by  Ivanov  et  al  (179) 
and  written  in  the  manuscript. 


(B.l) 
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B.2  Force  Expressions 


B.2.3  Hydrodynamic  Drag  Force  on  a  Single  Sphere  Trans¬ 
lating  Between  and  Parallel  to  Two  Parallel  Walls 

Detailed  multipole  solutions  for  the  drag  coefficient  for  a  single  solid  sphere  trans¬ 
lating  parallel  and  between  two  parallel  walls  was  given  by  Feuillcbois  et  al,  as 
cited  in  the  manuscript.  They  construct  an  interpolating  formulae  (eq.  29)  for 
this  drag  coefficient,  and  in  Figure  B.3  we  compare  the  results  of  this  correlation 
with  our  COMSOL  results  as  a  means  of  validating  our  numerical  simulations. 
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Figure  B.3:  Hydrodynamic  drag  force  on  a  single  sphere  translating  between  two 
parallel  walls  for  no  slip,  A  =  0  pm  as  a  function  of  the  sphere/wall  separation 
distance  d  (edge-to-edge),  comparing  the  Feuillebois  correlation  and  our  COMSOL 
calculation  for  two  values  of  the  channel  height  h  relative  to  the  sphere  radius  a. 
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